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ABSTRACT 



All Navy aircraft are required to accommodate, 
anthropometrically , ninety percent of the user population. 
Some designs have been criticized for their low accommo- 
dations but those accommodations have never been quantified. 
The purpose of this thesis was to quantify the accommoda- 
tion, by each type of operational Naval aircraft, of 
populations of Naval aviation personnel of 1964, 1969, and 
1975. The Computerized Accommodation Percentage Evaluation 
(CAPE) model was used to generate data points since only 
summary statistics were available for two of the popula- 
tions. Each subject of every population was checked against 
the requirements of the design specification, and against 
the limitations of each aircraft. All aircraft were found 
to accommodate more than ninety percent of the 1975 popu- 
lation. Time related changes in the pQpulations were noted 
and unexplained inconsistencies in the data were discovered. 
Possible sources of error were discussed and potential 
solutions proposed. 
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INTRODUCTION 



I . 

The Human Engineering Design Criteria for Military 
Systems, Equipment and Facilities [MIL-STD-1472B3 states 
the requirement for anthropometric accommodation of 
equipment users: 

Design and sizing shall insure accommodation, com- 
patibility, and maintainability by at least 90 
percent of the user population. Generally, design 
limits shall be based upon a range from the 5th 
percentile to the 95th percentile values for 
critical body dimensions. 

For the special case of Naval Aircraft cockpits, the 
accommodation criterion was enlarged in 1973 to include the 
central 95 percent of the user population CnAVAIR SD-24K, 
1973] . The anthropometric description of the population is 
based on data collected in 1964. Although a future survey, 
scheduled for 1981, will provide a more current description 
of the anthropometric features of Naval aviation personnel, 
aircraft now being designed for production and use well into 
the future are being designed to accommodate the 1964 popu- 
lation. The F-18, projected to become operational in 1982 
will be built to accommodate the pilots of 1964. Other 
designs, if they have advanced much beyond the concept 
stage, are committed to the same aircrew accommodation 
standards . 
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The anthropometric characteristics of the population 
are changing but the design specification is not. That the 
population is changing is reflected in efforts by the 
National Aeronautics and Space Administration and by the 
United States Air Force to predict the typical pilot size 
of the future. A 1978 study conducted by NASA concluded 
that stature (height) increases among pilot and potential 
astronaut males could be expected to amount to about 8 mm 
(1/3 inch) per decade [NASA 1024 Vol. I, 19781. The 
U.S. Air Force predicted an increase in height of 0.418 
inch per decade Roebuck, Kroemer, and Thompson, 1975 . 

These changes, since 1964 and into the future, affect 
pilot accommodation in operational Naval aircraft, but the 
exact nature of the effect is not known. 

In 1977 the Navy implemented the Anthropometric 
Compatibility Assignment Program COPNAVINST 3710.363, 3- 
system aimed at preventing the assignment of Naval Aviators 
and Naval Flight Officers to aircraft which are incompatible 
with their anthropometric features. Under the provisions of 
this program, pilots' and NFOs ' anthropometric measurements 
are taken when they enter flight training CbUMEDINST 3710.1, 
1977]. Values for sitting height, functional reach, 
buttock-knee length, and leg length (for descriptions of 
anthropometric terms, see Appendix A) are encoded and 
become part of their permanent records. For example, a 
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man coded 



6 4 2 3 

would have a sitting height measurement between 38.5 and 
38.9 inches, a functional reach of 30.5 to 30.9 inches, 
buttock-knee length of 23.0 to 23.9 inches and a leg length 
of 40.0 to 42.9 inches. 

The cockpits of each operational aircraft type are 
measured to determine the limits of pilot size compatible 
with the cockpit CnAVAIRINST 3710.9, 19793. The maximiom 
and minimum acceptable dimensions for sitting height, 
functional reach, buttock-knee length, and leg length are 
determined for each cockpit. These maxima and minima are 
then converted to four coded descriptors, which signify 
dimensions of exclusion for pilots. For example, an A-4E 
is coded; 

089 9 789 0 

As with the personal anthropometric codes, the first of the 
four groups of numbers corresponds to sitting height. The 
second group corresponds to functional reach; the third to 
buttock-knee length, and the fourth to leg length. The 
code 089 9 789 0 indicates that a pilot with a sitting 
height code 0, 8, or 9 would not be eligible for assignment 
to an A-4E. Similarly any pilot with a functional reach 
code 9, buttock-knee length code 7, 8, or 9, or a leg 
length code 0 would be excluded. 

The Chief of Naval Aviation Training (CNATRA) is 
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required to utilize these codes early in each aviation 
student's training to exclude the student from training 
sequences leading to aircraft with which he is not com- 
patible CCNATRAINST 13520.1, 1980H . The Navy Military 
Personnel Command uses the anthropometric codes in the 
process of assigning pilots and NFOs to aircraft squadrons. 

Currently, a similar coding scheme is under consider- 
ation for use internationally among free world nations 
CClauser, 1980] . The need for such a program was indi- 
cated by difficulties encountered by some foreign countries 
operating aircraft acquired from the United States and 
other foreign sources. 

In Figure 1, it is apparent that the height of U.S. 

Air Force flight personnel is generally greater than that 
of other countries. The fifth percentile in U.S. height 
is greater than the 50th percentile of the samples from 
Japan, Thailand, and Vietnam. Figure 2 shows a similar 
relationship for sitting height. 
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Figure 1 

Heights of Seven Military Populations 
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Figure 2 

Sitting Heights of Seven Military Populations 
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In this case the 5th percentile is larger than the 50th 
percentile of Thailand and Vietnam. Aircraft designed to 
accommodate people as small as the 5th percentile of the U.S. 
population can be expected to exclude about half of certain 
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foreign populations and as much as 75 percent of the 
Vietnam aviation personnel [Kennedy, 1973]. 

Although design criteria require accommodation of the 
central ninety percent of the Naval aviation population, 
the existing population is imperfectly described by the 1964 
sample. Changes in the selection process, trends in the 
United States general population were cited as factors con- 
tributing to the differences between the 1964 and the 1969 
populations [Moroney, Kennedy, Gifford, and Provost, 1971]. 
In addition, designing to accommodate the fifth and ninety- 
fifth percentiles in each variable does not yield an 
accommodation of ninety percent of the total population. 
Since exclusion results from any one measurement being out- 
side limits, accommodation is as much a function of the 
correlation between measurements as it is a function of 
each measurement. The correlation between measurements is 
appreciably less than 1.0 and accommodation is less than 
that indicated by any univariate estimate. When twelve 
measurements were considered, the exclusion of subjects who 
were outside the 5th and 95th percentile range in any 
measurement yielded a cumulative exclusion of 48 percent 
[Moroney and Smith, 1972]. As the existing population is 
altered further by expanding the acceptable range to 
include a greater percentage of the female population, user 
accommodation will be adversely affected. For example, a 
generalized cockpit built to male anthropometric standards 
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excluded 22 percent of the male user population. The same 
cockpit excluded 90 percent of the females in the U.S. Air 
Force CKetcham-Weidl and Bittner, 19761. 

Before considering future changes in the population, 
the current anthropometric exclusion situation should be 
determined. However, the proportion of the Naval aviation 
population excluded by the existing design specification 
has not been quantitatively determined. Similarly the 
proportion excluded by each specific aircraft design is 
unknown. Some aircraft can be expected to exclude more 
than the ten percent designed out by the specification. 

Other aircraft will exclude less than the maximum allowable 
ten percent by virtue of design allowance or perhaps by 
chance . 

The methods employed in measuring aircrew personnel 
have been found to be rather unreliable CMoroney, et . al . , 
197lJ . These measurements, however, form the basis of the 
Anthropometric Compatibility Assignment Program and their 
accuracy is critical since they constitute a determining 
factor in an officer's career path starting point. Addition- 
ally, these same data will be used to help identify (or 
eliminate) subsequent aviation duty assignments. Consider- 
able attention has been placed on determining the proper 
use of the data but the actual collection of the data 
appears to be an area where reliability could be improved. 
Currently, there is no systematic method of checking 
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anthropometric measurements. Each student is measured once 
without any attempt at remeasurement or checking of the 
accuracy of the measurements. If a measurement is recorded 
in error, the error becomes part of the student's record, 
undetected until some independent event brings the error to 
light. Errors in measurement of as much as three inches 
have been found Co' Leary, 1980] subsequent to the erroneous 
measurement being used to determine anthropometric suit- 
ability and assignment of a student to an aircraft type. 

Some students have completed flight training before anthro- 
pometric measurement errors were detected, resulting in 
reassignment and expensive retraining. 

The purpose of this thesis is to quantify the amount 
of incompatibility between aircraft and aircrewmen and to 
propose some alternatives to improve the reliability of the 
personal anthropometric measurements now being taken. By 
examining samples drawn from three Naval aviation popula- 
tions, each of a different year, proportions of those 
populations excluded from Naval aircraft will be determined. 
The exclusions resulting from changes in the population will 
be quantified and the specific proportions excluded from 
each aircraft type, as a function of sample source, will be 
calculated. Those results will lead to a discussion of some 
inconsistent patterns in the data and recommended changes in 
the data collection procedure now used by the Navy. 
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Chapter 2, DATA SOURCES, describes the three data sets 
analyzed. For each of the three samples: the subjects were 

described, previous statistical analyses discussed, and 
methods of data collection compared. Chapter 3, DATA 
PREPARATION, describes the preparation of each data set 
to facilitate analysis. Preparations consisted primarily 
of the use of a segment of the Computerized Accommodation 
Percentage Evaluation (CAPE) model to generate data points 
from a statistical description of the population. Chapter 
4, ANALYSIS OF DATA, describes the screening of each data 
set to determine the proportion of the sample excluded. 
Screening involves the elimination of any subject with an 
anthropometric measurement outside prescribed limits. The 
proportions excluded from each sample are then compared. 
Observations on the apparent shortcomings in data collection 
methods, errors in measurements, and unexplained inconsis- 
tencies were included. Chapter 5, SUMMARY AND CONCLUSIONS, 
includes inferences concerning interpretations of anthropo- 
metric data and recommendations in the area of data 
collection methods for use by the Navy. 
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I I . DATA SOURCES 



The data to be analyzed consisted of samples from three 
populations. Samples of Naval aviation personnel were taken 
in 1964, 1969, and 1975. Differences between populations 
other than the year of the sample exist and they will be 
discussed below. 

A. THE 1964 SURVEY DATA 

In 1964, a survey was conducted to determine the anthro- 
pometric characteristics of personnel flying in Naval 
aircraft. The survey was conducted in a manner to yield 
results sufficiently precise to be utilized subsequently for 
design of aircraft. The Bureau of Medicine and Surgery, 
tasked to collect data for use by the Bureau of Weapons, 
sent a trained anthropometric measuring team to ten different 
Naval and Marine Corps air stations to take measurements. 

Included in the sample were pilots. Naval Aviation 
Observers, bombardier-navigators, radar observers, flight 
surgeons who were designated Naval Aviators, and enlisted 
personnel who were permanently designated aeronautical 
personnel and who were still eligible for assignment to 
subsequent duty involving flying. Naval Aviators included 
both U.S. Navy and the U.S. Marine Corps flyers. 

The subjects in this survey comprised approximately ten 
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percent of the Naval aviation population on active duty in 
1964. The sample included personnel of various ages, ranks, 
ethnic backgrounds and sources (Naval Academy, NROTC, etc.). 
In that respect the sample undoubtedly provided an accurate 
representation of the Naval aviation population. The com- 
position of the sample is somewhat different from the later 
samples in the respect that enlisted men and older men were 
surveyed in 1964 but not in the later samples. The relia- 
bility of the measurements received special attention and 
was enhanced by the repeat measurement of randomly selected 
subjects, insuring rapid detection of measurement errors and 
encouraging the data collectors to be careful in their task. 
Each measurement was taken by personnel specifically trained 
for the job, using a technique of measurement recognized as 
the most accurate available. 

Published results of the survey include means and 
standard deviations, as well as percentile listings for each 
variable COifford, Provost, and Lazo, 19653. Subsequent 
analysis of the data led to publication of correlations 
between pairs of variables CMoroney, et . al . , 1971]. 

The results of this survey were adopted by the Navy as 
the basis for aircraft design specifications. MIL-STD-1333A , 
Military Standard Aircrew Station Geometry for Military 
Aircraft, includes NAEC-ACEL-553 as part of the standard, 
thus specifying the 1964 data as the guideline for aircraft 
being built today (NAEC-ACEL-553 reported the results of the 
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1964 survey). In that sense, these results are recognized 
as the standard and they have become the yardstick against 
which other data can be compared. 

This data set has two shortcomings with respect to this 
thesis, however. No measurements of leg length (or buttock- 
leg length) were taken. The components of leg length 
(buttock-knee length, knee height sitting and popliteal 
height sitting) were recorded but leg length had not yet 
been recognized as a critical selection dimension. Also, 
the results of data analysis are available, but the raw 
data are not. The importance of these two limitations will 
become apparent when procedures are discussed. 

B. THE 1969 DATA 

Subsequent to the introduction of an integrated anthro- 
pometric measuring device in 1964, the Navy required that 
the device be used to measure all aviation training candi- 
dates during their aviation training entrance physical 
examination, which was administered at the School of Aviation 
Medicine (now the Naval Aerospace Medical Institute). The 
collection of anthropometric data was thus institutionalized 
and became a routine process. Eight measurements were taken 
on each entering Student Naval Aviator and Student Naval 
Flight Officer. These measurements were: 

weight trunk height 

height buttock-knee length 

sitting height buttock-leg length 

shoulder width functional reach 
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Between January 1966 and August 1969, 6,534 students 
were measured. Their measurements comprise the data set 
herein called the 1969 data. The description of this data 
set was published in 1971 CMoroney, et . al . , 19711 and 
includes the mean and standard deviation of six measurements 
(shoulder height and functional reach were deleted from the 
analysis), percentile listings for those measurements and 
correlations between variables. 

No information on functional reach or shoulder height 
was published because of the discovery of measurement tech- 
nique errors resulting in measurements uncorrectible by 
any means short of repeat measurement. This deviation from 
ideal procedure should not have been altogether surprising 
since measurements were taken by personnel whose level of 
training did not match that of the 1964 data collectors. 

In addition, measurements were taken and recorded with no 
subsequent checking, leaving procedural errors undetected 
and biasing the measurements of many subjects. 

The size of the 1969 data set is quite large, consisting 
of virtually all of the student aviation input for over three 
years. As with the 1964 data set, no raw data is available, 
only analytical results. The absence of functional reach 
from the data analysis will present an obstacle to be 
discussed in Chapter 3. 
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C. THE 1975 DATA 



Following the publication of the results of the 1969 
data analysis, an improved integrated anthropometric 
measuring device was installed at NAPvII and it was subse- 
quently used to measure incoming students. That device is 
still in use today. It is operated by NAMI staff Hospital 
Corpsmen and it was used to measure the same eight dimen- 
sions as was its predecessor. A sample of 968 students 
who entered training in 1975 comprise the subjects of the 
1975 data. The raw data is available and from an examina- 
tion of the data points thirty-two subjects were found to 
be deficient in some respect and did not meet the entry 
requirements delineated by the Bureau of Medicine and Surgery 
LMANMED, 1978]. Most of these ineligibilities were based on 
excessive weight. However, all subjects were screened for 
minimum and maximum values of weight, height, sitting height, 
and buttock-leg length, as well as for height-weight rela- 
tionship. The remaining 936 subjects comprise the acceptable 
sample . 

The 1975 data were collected in a manner similar to that 
used for the 1969 data set; a routine execution of estab- 
lished procedures by Navy enlisted personnel. The 1975 
data also manifest data collection deviations similar to 
those of the 1969 data. The number of subjects is smaller 
than either of the other samples but still large enough to 
lend itself to statistical analysis. 
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III. DATA PREPARATION 



The goal of the analysis was to determine the proportions 
of the 1964, 1969, and 1975 samples excluded by the four 
anthropometric limitations (sitting height, functional 
reach, buttock-knee length, and leg length) required by 
OPNAVINST 3710. 36A. This goal was achieved by screening 
each subject four times, once for each of his measurements, 
and eliminating every subject whose measurements were not 
all within limits. 

This procedure required knowledge of the four individual 
measurements for each subject. In the 1975 data set, the 
four measurements (and more) were readily available. For 
the earlier data sets, only means, standard deviations, and 
correlations were available and no individual measurements 
were known. It was necessary, therefor, to generate data 
to simulate those measurements. 

The Computerized Accommodated Percentage Evaluation 
(CAPE) model C Bittner, 19751 was used to generate the 
required data points. This model required as entering 
arguments the mean and standard deviation for each of the 
four variables, as well as a matrix of correlations between 
variables. The model used a Monte Carlo process to generate 
data points. A program listing is included in Appendix C 
(Appendix D presents the same program modified for use with 
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the IBM 360-67 computer and other machines using Fortran IV). 
The CAPE model has been incorporated in the Computerized 
Accommodated Reach Model and is available at the Naval 
Aviation Development Center, Warminster, Pennsylvania. 

Unfortunately, some of the values required for generating 
the data points were missing from the survey reports. 
Specifically, the 1964 data included no information on leg 
length and the 1969 data omitted functional reach. It was 
necessary to find values in similar populations and substitute 
them for those missing. Taking values from other populations 
required the assumption that concurrent populations, similar 
in cultural, racial, and demographic composition, and having 
passed similar anthropometric entrance requirements, will 
have comparable anthropometric features. That assumption 
was made and surrogate values for the missing means and stan- 
dard deviations were taken from populations which matched 
the descriptions of the populations under analysis (1964 and 
1975), both subjectively and numerically. Inter-variable 
correlations were determined on much the same basis: substi- 

tuting correlations from similar populations. Ketcham-Weidl 
and Bittner Ll9763 had substituted correlations from a data 
set of male subjects into a set of female data and exper- 
ienced little loss of precision. 

For the 1964 data, the mean and standard deviation for 
leg length were assumed to be 43.82 and 2.01 inches, respec- 
tively. Those values were established in a 1968 survey of 
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U.S. Army helicopter students CSchane, Littell and Moultree, 
1969]. Of the information available, the population of that 
survey, better than any other, approximated the 1964 Naval 
aviation population with respect to the factors which Moroney 
Cl97ll reported as influencing differences between the 1964 
and 1969 Navy data sets. The Sitting height-Leg length cor- 
relation maintained the same value (0.45) for the 1969 and 
1975 data sets so the same value was used for 1964. Func- 
tional reach-Leg length and Buttock knee-Leg length 
correlations (0.639 and 0.860) were extracted from a 1970 
survey of U.S. Army aviators CChurchill, McConville, Laubach, 
and White, 1971], again because of the similarities in 
descriptions of the 1970 Army and 1964 Navy populations. 

For the 1969 data. Functional reach values for USAF 
flight personnel reported by Churchill Cl97l] were used for 
the mean and standard deviation (31.24 and 1.62 inches). 

The sitting height-functional reach and the buttock knee- 
functional reach correlations were assumed to be the same as 
those found in the 1975 data, 0.44 and 0.63, respectively. 
These values were greater than those of the 1964 data and 
their use, rather than the 1964 values, would influence 
error in the conservative direction, that is, in the direc- 
tion of lesser difference between univariate and multi- 
variate exclusion. The leg length-functional reach corre- 
lation was taken from Churchill's 1970 survey and is the 
same value assumed for the 1964 data. 
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Table I 



CAPE INPUTS AND OUTPUTS 
FOR THE 1964 DATA 

Means, Standard Deviations and Ranges 



Variable 


Observed 
~ Data 


Generated 
" Data 


Sitting Height 
mean 
std.dev. 
range 


36.28 

1.25 

32.19-41.62 


36.29 

1.26 

31.81-40.25 


Functional Reach 
mean 
std.dev. 
range 


31.51 

1.42 

27.26-36.31 


31.47 

1.40 

26.35-35.38 


Buttock-Knee Length 
mean 
std.dev. 
range 


24.09 

1.00 

20.73-27.81 


24.10 

9.94 

21.06-27.11 


Leg Length^ 
mean 
std.dev . 
range 


43.82 

2.01 


43.82 

1.99 

37.97-49.87 



^Schane , W.P. , Littell, D.E.', and Moultree, C.G. , 
Selected Anthropometric Measurements of 1,640 U.S. Army 
Warrant Officer Candidate Flight Trainees , U.S. Army 
Aeromedical Research Laboratory Report 69-2, as reported 
in : Anthropometric Source Book Volume II: A Handbook of 

Anthropometric Data, NASA Reference Publication 1024, 
1978, p. 156. 
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Table II 



CAPE INPUTS AND OUTPUTS 
FOR THE 1964 DATA 

Correlations 





Sitting 

Height 


Functional 

Reach 


Buttock- 
Knee Length 


Leg 

Length 




Observed 


L Data 






Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


1.0 


.376 

1.0 


.382 

.586 

1.0 


.45^ 

.639t> 

.86b 

1.0 


Generated Data 


Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


1.0 


.410 

1.0 


.374 

.587 

1.0 


.446 
.650 
. 867 
1.0 



^Same as observed value for 1969 and 1975 samples. 

^Churchill, E. , McConville, J.T. , Laubach, L. , and 
White, R.M. , Anthropometry of U.S. Army Aviators-1970 , U.S. 
Army Natick Laboratories Technical Report 72-52-CE, December 
1971, p. 278. 
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Table III 



CAPE INPUTS AND OUTPUTS 
FOR THE 1969 DATA 

Means, Standard Deviations, and Ranges 



Variable 


Observed 

Data 


Generated 

Data 


Sitting Height 
mean 
std . dev. 
range 


36.76 

1.21 

31.70-41.0 


36.77 

1.22 

32.43-40.61 


Functional Reach^ 
mean 
std . dev . 
range 


31.24 

1.62 


31.20 

1.60 

20.02-35.73 


Buttock-Knee Length 
mean 
std . dev. 
range 


24.54 

1.26 

20.40-29.90 


24.55 

1.25 

20.67-28.31 


Leg Length 
mean 
std . dev. 
range 


43.86 

2.08 

32.50-50.70 


43.86 

2.05 

38.14-50.02 



^Churchill, E. , McConville, J.T. , Laubach, L. , and 
White, R.M. , Anthropometry of U.S. Army Aviators-1970 , U.S. 
Army Natick Laboratories Technical Report 72-52-CE, December 
1971, p. 91. 
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Table IV 



CAPE INPUTS AND OUTPUTS 
FOR THE 1969 DATA 

Correlations 





Sitting 

Height 


Functional 

Reach 


Buttock- 

Knee-Length 


Leg 

Length 




Observed 


[ Data 






Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


1.0 


, 44 a 

1.0 


.40 

.63^ 

1.0 


.45 

.639^^ 

.86 

1.0 


Generated Data 


Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


1.0 


.473 

1.0 


.393 

.629 

1.0 


.449 

.652 

.801 

1.0 



^Same as observed value in 1975. 

^Churchill, E. , McConville, J.T. , Laubach, L. , and 
White, R.M., Anthropometry of U.S. Army Aviators-1970 , U.S. 
Army Natick Laboratories Technical Report 72-52-CE, December 
1971, p. 278. 
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Table V 



DATA SET COMPARISON 



Means, Standard Deviations, and Ranges 



Variable 


1964 


1969 


1975 


Sitting Height 


mean 


36.28 


36.76 


36.61 


std . dev . 


1.25 


1.21 


1.15 


range 


32.19-41.62 


32.43-40.61 


32.90-40.00 


Functional Reach 


mean 


31.51 


31.24 


30.89 


std . dev . 


1.42 


1.62 


1.16 


range 


27.26-36.31 


26.02-35.73 


27.10-34.60 


Buttock-Knee Length 


mean 


24.09 


24.54 


24.05 


std. dev. 


1.00 


1.26 


1.07 


range 


20.73-27.81 


20.67-28.31 


20.70-26.80 


Leg Length 


mean 


43.82 


43.86 


42.89 


std . dev . 


2.01 


2.08 


1.89 


range 


37.97-49.87 


38.15-50.02 


37.10-48.20 
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Table VI 



DATA SET COMPARISON 
Correlations 



1964 





Sitting 

Height 


Functional 

Reach 


Buttock- 
Knee Length 


Leg 

Length 


Sitting Height 


1.0 


.376 


.382 


.45 


Functional Reach 




1.0 


.586 


.639 


Buttock-Knee Length 






1.0 


.86 


Leg Length 








1.0 



1969 



Sitting Height 1.0 


.44 


.40 


.45 


Functional Reach 


1.0 


.63 


.639 


Buttock-Knee Length 




1.0 


.86 


Leg Length 






1.0 



1975 



Sitting Height 1.0 


.44 


.40 


.45 


Functional Reach 


1.0 


.63 


.61 


Buttock-Knee Length 




1.0 


.79 


Leg Length 






1.0 
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When the assumed values described above were included 



in the data analysis, the requirements for CAPE input were 
complete. The inputs to CAPE are listed in Tables I 
through IV under Observed Data, One thousand subjects were 
generated for the 1964 and for the 1969 data sets. Together 
with the 1975 data in its original, as observed form, the 
result was three similar data sets with sample sizes of 
1000, 1000, and 936. The validity of this model's perfor- 
mance has been demonstrated with a correlation between 
empirical and generated results of 0,997 to 0.999 CBittner, 
19743, For the 1964 and 1969 data sets, the CAPE generated 
output was analyzed to insure that the means, standard 
deviations, and correlations of the generated data were 
actually those sought. The values used as inputs to CAPE 
and those resulting from analysis of the generated data are 
shown in Tables I through IV. Tables V and VI compare the 
summary statistics for all three data sets. 
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IV. ANALYSIS OF DATA 



As a result of the data generating procedure, all three 
data sets included a value for each of the four relevant 
measurements for each subject. The samples were then 
examined more closely to determine which subjects in each 
sample would be excluded, either by design constraints (5th 
to 95th percentiles or 3rd to 98th percentiles) or by air- 
craft constraints specified by NAVAIRINST 3710.9. 

A . PERCENTAGES EXCLUDED 

Fortran program SCREEN was used to count the subjects 
who were too large or too small to satisfy anthropometric 
requirements. For each of the three samples, identical 
screening operations were conducted. Each subject was 
examined and excluded if any one or more of the four 
measurements under consideration was outside prescribed 
limits. Each subject's anthropometric features were tested 
against each cutoff limit. Thus an individual might be 
excluded on both functional height and sitting height but 
not on buttock-knee length and leg length. For each elimin- 
ation encountered, the subject was counted as being 
eliminated and each variable causing an elimination was 
tallied. At the end of each sample, the resultant totals 
were converted to percentages of the total sample, yielding 
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a percentage eliminated by each of the four variables and 
a total exclusion resulting from the net effect of the four 
variables' eliminations. 

The prescribed limits mentioned above were set equal to 
the percentile values in the design specification (5th and 
95th or 3rd and 98th percentiles). Then the limits were 
set, successively, to the anthropometric limits of each 
aircraft for which anthropometric coding had been estab- 
lished CnAVAIRINST 3710.9]. 

1 . Design Constraints 

To screen the samples by design specification, the 
subjects were screened to eliminate those with any of the 
four measurements less than the 5th percentile or above the 
95th percentile of the 1964 data, the design standard CmIL- 
STD-1472B]. Then the screening cutoff points were changed 
to correspond to the 3rd and 98th percentiles CnAVAIR SD-24K, 
1973], and the screening process was repeated. 

It is interesting to note at this point that percen- 
tile values for leg length of the 1964 data do not exist 
since that measurement was not taken in the 1964 survey. 
Still, the Navy considers leg length to be a critical dimen- 
sion in that it is one of the measurements included in the 
Anthropometric Compatibility Assignment Program. Critical 
dimensions are also the basis of the design specification 
since '...design limits shall be based upon a range from the 
5th percentile to the 95th percentile values for critical 
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body dimensions' CMIL-STD-1472B3 . Since the percentile 
values were not available, values used for leg length 
percentiles were determined from the CAPE generated sample 
of 1000 subjects, described in Chapter III, DATA PREPARATION. 

The design specification percentile screens resulted 
in the maximum exclusion percentage to be expected if air- 
craft were designed to accommodate only the specified 
percentiles, and no more, for each dimension. This result 
amounts to a worst case exclusion and could occur only if 
a design exactly met the bare minimum accommodation in each 
variable, a very unlikely event, and one which would still 
not meet the requirement to accommodate 90 percent of the 
user population as specified by MIL-STD-1472B. The percen- 
tile screen results, shown in Table VII, show that cumulative 
exclusion can greatly exceed the maximum single variable 
elimination. For example, the 1969 data screened on the 
5th to 95th percentiles shows a total exclusion of 40.4 
percent but the highest single variable elimination is the 
21 percent excluded by buttock-knee length. 
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Table VII 



Design Specification Exclusions for 
the 5-95 Percentiles and the 3-98 Percentiles 



Variable 


1964 


1969 


1975 




5th to 


95th Percentile 




Percentage Excluded 


Sitting Height 


11.1 


11.5 


7.7 


Functional Reach 


8.7 


14.6 


8.4 


Buttock-Knee Length 


9.1 


21.0 


13.2 


Leg Length 


9.8 


11.3 


11.8 


Cumulative 


26.4 


40.4 


26.7 




3rd to 


98th Percentile 




Percentage Excluded 


Sitting Height 


5.0 


4.7 


2.7 


Functional Reach 


4.8 


9.4 


5.9 


Buttock-Knee Length 


5.4 


11.6 


6.4 


Leg Length 


4.9 


4.7 


6.9 


Cumulative 


14.3 


22.7 


14.3 
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The difference between cumulative exclusion and the 



largest univariate exclusion figure, which could be taken 
falsely as an upper bound, is due to the non-unity corre- 
lation between variables. Although some anthropometric 
models have described all anthropometric dimensions as a 
specific multiple of height I^Roozbazar, 19793 with the 
implicit assumption of intervariable correlation equal to 
1.00, those models fail to recognize the variable propor- 
tions of people and would lead to the acceptance of the 
largest univariate exclusion as the net cumulative exclu- 
sion. The assumption of a correlation of 1.00 amounts to 
treating the population as if the man with 95th percentile 
height must also have 95th percentile leg length, 95th 
percentile shoulder width, and 95th percentile in all 
measurements, i . e . , 'the 95th percentile man.' If that were 
the case, the percentile screens would have eliminated the 
same subjects on all variables, as well as on the cumulative 
exclusion. If, at the other extreme, intervariable corre- 
lation were zero, the single variable eliminations would be 
independent and their arithmetic sum would appear as the 
cumulative exclusion. 

In reality, some subjects with leg lengths shorter 
than the 5th percentile also have sitting heights less than 
the 5th percentile. Others have short legs but sitting 
heights nearer the mean. The evidence of variable propor- 
tions along with variable sizes is found in the cumulative 
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exclusion values. Generally, the cumulative exclusion is 
less than the sum of the single variable eliminations but is 
greater than any of the single value eliminations. Only 
rounding error could cause the cumulative exclusion to 
exceed the siim of the single variable eliminations. The 
cumulative exclusion is generally as small as the largest 
single variable when only one variable resulted in any 
elimination . 

2 . Aircraft Constraints 

The specific aircraft screens were accomplished with 
the same technique as used on the design specification 
screens. The cut points for each variable were the values 
corresponding to each aircraft's anthropometric restrictions 
specified by NAVAIRINST 3710.9. Then, for each aircraft, 
the three samples were examined to eliminate subjects whose 
anthropometric measurements were large enough, or small 
enough, to be excluded by the aircraft's restrictions. The 
number of subjects eliminated by each measurement and the 
cumulative exclusions were totaled and converted to a percen- 
tage of the sample. Different aircraft with the same 
anthropometric restrictions were consolidated into eighteen 
groups and the screens were executed, resulting in the 1964, 
1969, and 1975 samples being matched to all aircraft types. 
The percentage exclusions are shown in Table VIII, which 
lists aircraft types in alphabetical order. 
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Table VIII 



Exclusion 


by Aircraft Type 


Aircraft Code Variable 


Exclusion Percentage 
1964 1969 1975 



A-3 


— 


sit . ht . 


— 


— 


— 




9 


fnct . rch . 


0.6 


2.6 


0.3 




— 


b-knee Inth. 


— 


— 


— 




0 


leg Inth. 
cumulative 


0.1 

0.6 


0.0 

2.6 


0.1 

0.4 


A-4E 


089 


sit . ht . 


0.9 


1.7 


0.4 


A-4F 


9 


fnct . rch . 


0.6 


2.6 


0.3 




789 


b-knee Inth. 


1.0 


6.6 


0.6 




0 


leg length 
cumulative 


0.1 

2.5 


0.0 

10.3 


0.1 

1.5 


A-4M 


0 


sit . ht . 


0.2 


0.1 


0.0 




9 


fnct . rch . 


0.6 


2.6 


0.3 




789 


b-knee Inth. 


1.0 


6.6 


0.6 




0 


leg length 
ctimulat ive 


0.1 

1.8 


0.0 

9.3 


0.1 

1.1 


A-6A 


06789 


sit . ht . 


5.0 


9.6 


4.0 


A-6B 


— 


fnct . rch. 


— 


- 


— 


A-6E 


9 


b-knee Inth. 


0.0 


0.4 


0.0 






leg length 
cumulative 


5.0 


9.9 


4.0 


A-7A 


089 


sit . ht . 


0.9 


1.7 


0.4 


A-7C 


6789 


fnct . rch . 


7.4 


13.6 


8.8 


A-7E 


9 


b-knee Inth. 


0.0 


0.4 


0.0 




0 


leg length 
cumulative 


0.1 

8.2 


0.0 

15.6 


0.1 

9.2 


AV-8A 


09 


sit . ht . 


0.4 


0.8 


0.1 




9 


fnct . rch. 


0.6 


2.6 


0.3 




9 


b-knee Inth. 


0.0 


0.4 


0.0 






leg length 
cumulative 


1.0 


3.8 


0.4 


C-118 




sit . ht . 


__ 


_ 


_ 


C-131 


— 


fnct . rch . 


- 


- 


- 






b-knee Inth. 


— 


- 


— 




0 


leg length 
cumulative 


0.1 

0.1 


0.0 

0.7 


0.1 

1.6 
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Table VIII (Continued) 
Exclusion by Aircraft Type 



Exclusion Percentage 

Aircraft Code Variable 1964 1969 1975 



CH-53 


01 


sit • ht . 
fnct . rch. 
b-knee Inth. 
leg length 
cumulative 


EA-3 


9 

0 


sit . ht . 
fnct . rch. 
b-knee Inth. 
leg length 
cumulative 


EA-6B 


06789 

9 


sit . ht . 
fnct . rch. 
b-knee Inth. 
leg length 
cumulative 


EA-6B 


06789 


sit . ht . 


(LF) 


9 

9 


fnct . rch . 
b-knee Inth. 
leg length 
cumulative 


F-4B(F) 


0789 


sit . ht . 


F-4J(F) 


89 


fnct . rch . 


F-4N(F) 


9 


b-knee Inth. 


F-4S(F) 


0 


leg length 
cumulative 


F-4B(R) 


89 


sit . ht . 


F-4J(R) 


- 


fnct . rch. 


F-4N(R) 


9 


b-knee Inth. 


F-4S(R) 


— 


leg length 
cumulative 



0.7 


0.7 


1.6 


0.7 


0.7 


1.6 


0.6 


2.6 


0.3 


0.1 


0.0 


0.1 


0.6 


2.6 


0.4 


5.0 


9.6 


4.0 


0.0 


1 o 


o 
1 o 


5.0 


9.9 


4.0 


5.0 


9.6 


4.0 


0.6 


2.6 


0.3 


0.0 


0.4 


0.0 


5.6 


12.5 


4.3 


2.0 


3.9 


1.7 


1.8 


4.3 


1.9 


0.0 


0.4 


0.0 


0.1 


0.0 


0.1 


3.8 


8.4 


3.6 


0.7 


1.6 


0.4 


O 1 

o 


0.4 


0.0 


0.7 


2.0 


0.4 
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Table VIII (Continued) 
Exclusion by Aircraft Type 



Exclusion Percentage 



Aircraft 


Code 


Variable 


1964 


1969 


1975 


F-8D 


0789 


sit . ht . 


2.0 


3.9 


1.7 


F-8J 


— 


fnct . rch. 


— 


- 


- 




9 


b-knee Inth. 


o 

o 


0,4 


o 

o 




- 


leg length 


- 


- 


- 






cumulative 


2.0 


4.2 


1.7 


F-14(F) 


0 


sit . ht . 


0.2 


0.1 


0.0 




9 


fnct . rch. 


0.6 


2.6 


0.3 




9 


b-knee Inth, 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


00 

o 


3.1 


O 

CO 


F-14(R) 





sit . ht . 












— 


fnct . rch. 


— 


— 


- 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


o 

o 


0.4 


o 

o 


H-2 





sit . ht . 











89 


fnct . rch. 


1.8 


4.3 


1.9 




— 


b-knee Inth. 


— 


— 


- 




0 


leg length 


0.1 


0.0 


0.1 






cumulative 


1.8 


4.3 


1.9 


OV-10 





sit . ht . 












9 


fnct . rch. 


0.6 


2,6 


0.3 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


0.6 


3.0 


O 

CO 


RA-5C 


09 


sit , ht . 


0.4 


00 

o 


0.1 


(F) 


— 


fnct . rch . 


- 


— 


- 




9 


b-knee Inth. 


0.0 


0.4 


0,0 




0 


leg length 


0.1 


0.0 


0.1 






cumulative 


0.5 


1.2 


0.2 


RF-4B 


0789 


sit . ht . 


2.0 


3.9 


1.7 


(F) 


89 


fnct . rch. 


' 1.8 


4.3 


1.9 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




0 


leg length 


0.1 


0.0 


0.0 






cumulative 


3.8 


8.4 


3.6 
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Table VIII (Continued) 
Exclusion by Aircraft Type 



Aircraft 


Code 


Variable 


Exclusion Percentage 
1964 1969 1975 


RF-4B 


89 


sit . ht . 


0.7 


1.6 


0.4 


(R) 


— 


f net . rch . 


— 


— 


- 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


0.7 


2.0 


0.4 


RF-8 


0789 


sit . ht . 


2.0 


3.9 


1.7 




— 


fnct . rch. 


— 


— 


— 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


2.0 


4.2 


1.7 


S-3 





sit . ht . 


_ 










— 


fnct . rch . 


— 


— 


— 




9 


b-knee Inth. 


0.0 


0.4 


0.0 




- 


leg length 


- 


- 


- 






cumulative 


0.0 


0.4 


0.0 


T-2C 


0 


sit . ht . 


0.2 


0.1 


0.0 


TA-4 


9 


fnct . rch. 


0.6 


2.6 


0.3 


TAV-8 


9 


b-knee Inth. 


0.0 


0.4 


0.0 




0 


leg length 


0.1 


0.0 


0.1 






cumulative 


0.8 


3.1 


0.4 


U-16 





sit . ht . 





_ 


_ 




— 


fnct . rch . 


- 


- 


— 




— 


b-knee Inth. 


— 


- 


- 




0 


leg length 


0.1 


0.0 


0.1 






cumulative 


0.1 


0.0 


0.1 



NOTE: 






(F) 


- Front 


seat 


(R) 


- Rear 


seat 


(LF) 


- Left 


front seat 
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B . RESULTS 



The design criterion requiring accommodation of the 5th 
and 95th percentiles results in exclusion of over 26 percent 
of the population. That result fails to meet the goal of 
90 percent accommodation. Designing to the 3rd and 98th 
percentiles, instead of the 5th and 95th, reduced exclusion 
to 14 percent. That result still does not meet the 90 per- 
cent accommodation goal. 

No aircraft type's coding eliminated as much as ten per- 
cent of the 1964 sample. The A-7 exclusion was 8.2 percent 
and was the worst case. The accommodation by all aircraft 
of more than 90 percent would seem to indicate that the 
cockpit designs included a margin of accommodation allowance 
in addition to the minimum required. Such a safety factor 
in design is not unexpected but is still technically an 
unrequired bonus. 

The exclusion percentage was not constant, but changed 
from one sample to the next. The greatest exclusions were 
found in the 1969 data while the 1964 and 1975 cumulative 
exclusions were about the same. Although cumulative exclu- 
sions were nearly equal, the contributing factors had 
varying proportions. Sitting height eliminations decreased 
and eliminations due to buttock-knee length and leg length 
increased. 

The 1969 data exhibited other features of incongruity 
which detracted from the credibility of the three data sets 
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Figure 3 



Comparison of Selected Variables 
and Standard Deviations 



Variable Inches 



Year 



1964 1969 1975 



Height 



73- 

72- 



70- 

69- 

68 - 

67- 



Leg 

Length 



46- 

45- 

44- 

43- 

42- 

41- 



No 

Data 



Sitting 

Height 



38- 

37- 

36- 

35- 



1 



1 



1 



Buttock- 
Knee Length 



27- 

26- 

25- 

24- 

23* 

22 - 






i 



1 



NOTE: o - mean 



T 



- one standard deviation from mean 
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taken as a whole. These effects reflected the presence of 
unidentified factors, not entirely explained by the facts 
available. While height increased steadily from 1964 to 
1975, its components did not. As seen in Figure 3, which 
shows mean values and one standard deviation on either side 
of the mean, mean height increased from 69.94 to 70.15 
inches, a rate of 0.427 inches per decade. That rate con- 
firms the USAF prediction of 0.418 inch per decade C Roebuck, 
et . al . 3 . However, the increased height does not appear as 
an increase in the other long bone measurements which are 
components of height. From 1969 to 1975, both mean sitting 
height and mean leg length appear to have decreased. The 
increase in height should be traceable to an increase in 
sitting height, an increase in leg length, or both, but 
neither was noted. Explanations other than erroneous data 
are possible, but not supported by knowledge of the popula- 
tion. It could be hypothesized that an increasing 
concentration of fat in the buttocks could increase both 
sitting height and leg length measurements while not 
affecting height. This effect has not been observed. 

Another possible explanation, operator induced loss of 
precision, would be accompanied by increases in the variation. 
An examination of Table V and Figure 3 indicate that standard 
deviations remained relatively constant, indicating the 
absence of excessive random error by equipment operators. 
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The inclusion of a constant error in any of the measure- 
ments in one or more data sets could account for the 
apparent inconsistency. This kind of error should be 
precluded by the use of standard procedures and calibrated 
equipment. The 1969 data were collected from measurements 
taken with a device for which no calibration was provided. 
There was no calibration procedure or equipment associated 
with the integrated anthropometric measuring device until 
after the introduction of the improved device in 1971. 
Although a means of calibrating the improved anthropometric 
measuring device was provided after 1971, the calibration 
equipment is no longer available and the date of the last 
calibration is unknown. 

The largest contributor to the apparent error is the 
diminished leg length in 1975. The reduction of nearly an 
inch is difficult to believe. Leg length is also a diffi- 
cult measurement to take accurately CMoroney, 19803. The 
position of the subject being measured (sitting erect with 
his right leg maximally extended horizontally in front of 
him) is difficult for some individuals to attain. The 
inflexibility of some muscular legs prevents the leg from 
being straightened completely. Methods aimed at straighten- 
ing the leg can result in the hips rolling forward, causing 
the buttocks to move outward from the seat back. The 
measured leg length is then either too long (with the leg 
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straight and buttocks forward) or too short (with the leg 
slightly bent ) . 

The 1975 data also appeared to have been influenced by 
rounding. The frequency of whole inch and half inch measure- 
ments is greater than expected. The bias incurred in the 
mental rounding process is indeterminate and no effort was 
made to correct for it. The existence of the effect indi- 
cates a need for a procedural modification. 

The exclusion determined here were based on the consider- 
ation of only four variables. Although, of the available 
measurements, they were considered to be the most important 
with respect to accidents, injuries, and aircraft controlla- 
bility, CMoroney, 1980] they are not the only measurements 
relevant to accommodation. Administrative limitations set 
by NMPC established that no more than four variables would 
be used. The practice of collecting anthropometric data at 
the flight training entry point and the specific measurements 
taken there determined which measurements could be used. 

The inclusion of additional measurements, such as shoulder 
width, could only increase exclusions. The substitutions of 
other measurements could affect the exclusion results 
similarly . 
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V. SUMMARY AND CONCLUSIONS 



The Navy requires airplane cockpit designs to accom- 
modate ninety percent of the Naval aviation population. In 
support of that goal, designs are required to accommodate 
the 5th to 95th percentiles (or the 3rd to 98th percentiles) 
of critical body dimensions. The success or failure of 
aircraft designers in accommodating the population had not 
been quantified, heretofore. By examining anthropometric 
data from three populations (1964, 1969, and 1975), two of 
which were generated by the Computerized Accommodated Percen- 
tage Evaluation Model (CAPE), exclusions allowed by the 
design specification and exclusions associated with each 
aircraft type were determined. The numerical values corres- 
ponding to the percentiles named in the specification (5 to 
95 or 3 to 98) were compared to the measurements of each 
subject in each population to determine the proportion of 
subjects of each population having at least one measurement 
outside the limits of accommodation and therefor excluded. 

The results were the exclusions expected from designing to 
the specified percentiles. The aircraft exclusions were 
found by comparing aircraft anthropometric exclusion codes, 
assigned by NAVAIR, and the measurements corresponding to 
those codes, to the three populations. That process 
resulted in an exclusion percentage for each aircraft and 
each population. 
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It was found that : 



1. The design specification requiring accommodation of the 
5th to 95th percentiles excluded 26 percent of the specified 
population . 

2. Designing to the 3rd to 98th percentiles excluded 14 
percent of the population. 

3. Exclusions, both those allowed by the specification and 
those peculiar to each aircraft, varied over time. 

4. All aircraft accommodate more than the required ninety 
percent of the most recent population sample. 

5. The proportion of exclusion caused by a given variable 
changed over time. 

6. There are inconsistencies in the data. 

A. DESIGN SPECIFICATION 

The design specification should result in at least ninety 
percent of the user population being accommodated by each 
aircraft. Reference to the 5th and 95th percentiles and to 
the 3rd and 98th percentiles of critical body dimensions 
serves only to fog the issue. The imaginary '95th percentile 
man,' consisting of a 95th percentile measurement for all 
variables, is assumed from the requirement to design to his 
dimensions. The tacit assumption that anthropometric 
features can be assembled and generalized to form a 95th 
percentile man is false and misleading since it assumes a 
correlation of 1.00 between anthropometric features. The 
specification would be more effective if it were to require 
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accommodation of ninety percent of the population and remain 
silent on the accommodation of singular dimensions. Deletion 
of all reference to percentiles in the specification would be 
a productive change. However, anthropometric data should 
still be reported in percentiles for use in univariate 
applications . 

The Navy is currently enjoying the benefit of an unre- 
quired bonus in the area of user accommodation. Ninety 
percent accommodation was required, more than that is 
delivered, at least for the four variables examined, based 
on the 1975 data. Although there appears to be nothing to 
indicate that this bonus will not continue, there is like- 
wise nothing to guaranty that it shall. If a future aircraft 
design were successful in accommodating exactly the required 
ninety percent, the consequences would amount to exclusions 
greater than those being experienced on any current Navy 
aircraft and the design specification would still have been 
satisfied. Since there is feeling in the Navy that current 
user accommodation is insufficient and since the ninety 
percent accommodation requirement is being met, either the 
Navy must change the ninety percent requirement or adjust its 
attitude toward current accommodations. The answer may be 
found in partitioning the user population into groups, each 
of which currently flies in a particular type aircraft. 

Then specify the 'user population' to be one of those groups, 
rather than all flying personnel in the Navy, and require 
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accommodation of the entire group. This solution, while 
easy to conceive, would be difficult to administer. 

B. COCKPIT RESTRICTIONS 

The exclusions enumerated in Table VIII indicate that 
the Navy has a minimal problem, on the global level, with 
incompatible planes and personnel. The attention that the 
problem has received would hardly seem justified by the 
miniscule exclusions indicated by the analysis. The 
expense in terms of money, man-hours, and effort reflected 
in the production of instructions by OPNAV, BUMED, CNATRA, 
and NAVAIR, together with the resultant effort necessary to 
comply with their requirements must have been justified by 
a judgement that a problem existed but the numbers in Table 
VIII do not reflect the magnitude of the problem. A reason 
that more impressive numbers don't appear in the table 
might be in the aircraft codes' inaccuracy, or in their 
failure to account for influential variables. 

The accuracy of the aircraft codes is dependent on a 
judgement factor on the part of the people who measure air- 
craft cockpits. The anthropometric codes assigned to 
aircraft are based on cockpit measurements adjusted for 
clothing and equipment worn by the user. Those adjustments 
cannot allow for all variations. The codes assigned to 
aircraft can also be affected by forecasts of the consequences 
of the codes assigned. The prospect of excluding an unsuit- 
ably large portion of the user population could have 
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influenced aircraft exclusion code assignments in the direc- 
tion of minimal effect. 

The high cost of error does not appear in the table but 
undoubtedly accounts for the attention given to such unim- 
pressive exclusion figures. The cost of losing an airplane 
or an aircrew is considerable. That cost makes any erroneous 
assignment of personnel potentially disastrous and may 
account for the emphasis in the area of aircrew accommoda- 
tion . 

Additionally, anthropometric restrictions can affect an 
officer's career rather markedly by denying him assignment 
to the type aircraft he is motivated to fly. That fact 
makes the globally small exclusion problem a crisis on the 
individual level and explains why the attention level 
exceeds that expected from examining the exclusion figures. 

C. DATA INCONSISTENCIES 

The data indicate that the values associated with anthro- 
pometric features shifted over time. Although the exclusions 
by design specification in 1964 and those in 1975 were about 
the same (see Table VIII), the individual variables respon- 
sible for the exclusions were of varying influence. In 1964, 
sitting height caused more exclusions than it did in 1975, 
while buttock-knee length and leg length caused more in 1975 
than in 1964. The 1969 data show higher exclusion than 
either of the other two data sets. Sitting height, buttock- 
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knee length, and leg length all peaked in 1969, accounting 
for the corresponding peak in exclusions. 

The exclusion results varied erratically over time. 

There are inconsistent effects in the data and the lack of 
continuity cannot be explained with the information avail- 
able. The leg length measurement is the most likely source 
of the inconsistent effects and the discontinuance of leg 
length as a critical measurement might be beneficial if 
another measurement could be substituted. Buttock-knee 
length might be used as a substitute for leg length since 
the two variables are well correlated (r=.8). Alternatively, 
a measurement taken in a sitting position approximating a 
cockpit seating position could be substituted. Measuring 
from the heel to the seat back via the front edge of the 
seat would accomplish, more directly, the purpose presently 
being attempted by the leg length measurement. The problems 
encountered measuring leg length would be eliminated and the 
measurements of personnel and of aircraft cockpits would be 
more comparable and, consequently, more effective in matching 
personnel to airplanes. Such a measurement of 'functional 
leg throw' would be a special purpose dimension not included 
in previous surveys. Previous surveys, however, were not 
faced with the special problem now at hand and the cost of 
designing and implementing this new measurement could be 
recovered in the prevention of a single aircraft accident. 
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D. DATA LIMITATIONS 



Of all the anthropometric dimensions which might be con- 
sidered, the Navy collects data on eight anthropometric 
features (including weight) considered relevant to flight 
safety. Of those eight variables, only four are considered 
in the process of assigning personnel to aircraft. The 
limitation to four variables is not necessarily optimum with 
respect to flight safety since the administratively driven 
maximum of four variables was the effective constraint. The 
inclusion of additional variables, if they were found to be 
relevant, could only decrease accommodation further. The 
decrease could be dramatic if the additional variables were 
poorly correlated to those already considered. 

The problems arising from variables in combination (but 
not visible in single variable evaluation) have been recog- 
nized [CNATRAINST 13520.1, 19803 but no description of the 
interaction between variables is included in the Anthropo- 
metric Compatibility Assignment Program. For example, the 
minimum functional reach and leg length necessary to reach 
controls varies as the position of the seat changes. The 
position of the seat is determined largely by the sitting 
height of its occupant. An individual with a short sitting 
height must raise the seat high enough to allow vision over 
the nose of the aircraft. As the seat rises, it also moves 
aft slightly, causing the rudder pedals to be adjusted aft 
and also requiring a longer functional reach to enable the 
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man to reach the controls. Thus, sitting height determines 
in part the effective minima for leg length and functional 
reach values, making required leg length and functional 
reach distances functions of sitting height. The values used 
for anthropometric coding, however, are determined irres- 
pective of seat position changes. An additional 
anthropometric code should be established to identify those 
personnel who meet all the limits but whose combination of 
anthropometric features could create problems. 

The procedures used to take measurements in the 1969 and 
1975 data sets were deficient in some respects. Most notably, 
the failure to calibrate the integrated anthropometric 
measuring device was unfortunate. The apparent rounding of 
measurements also induced error. The calibration deficiency 
can be corrected by instituting a requirement to document 
periodic calibration and adjustments, as is required by the 
Preventive Maintenance System. The rounding error could be 
eliminated by converting the measuring device to provide 
digital displays of measurements. Better, but also more 
expensive, an automated measuring device could record 
measurements without requiring the operator to read and 
record numbers. Such a device would also prevent the oper- 
ator from rounding measurements or inserting any bias into 
the data. 

The necessity to draw on outside sources for data con- 
cerning leg length, a critical body dimension about which 
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there is no data in the design specification (1964 data), 
points out a basic deficiency in the specification. 

Designers do not consider leg length in the design of 
cockpits. Instead, they use buttock-knee length and pop- 
liteal height sitting, which are components of leg length. 

The assignment process would be more consistent if similar 
variables were used for exclusion codes. 

E. CONCLUSIONS 

The data collection system which supports the Anthropo- 
metric Compatibility Assignment Program consists of 
measurements taken on aircraft cockpits as well as the 
equipment and procedures used to measure aviation personnel. 
Many inconsistencies in personnel measurements have been 
noted even though those measurements are an essential ingre- 
dient in the successful operation of the program. The 
program will operate no more effectively than its weakest, 
or limiting, factor. That limiting factor appears to be the 
measurement of personnel. Inconsistencies in the data appear 
to be results of measurement bias. The appearance of larger 
personnel in 1969 and corresponding higher exclusions in 
that population are more likely the changes of measurement 
bias than any real effect in the population. 

Measurement error could be reduced by: 

1. Periodically calibrating the measurement equipment. 

2. Randomly remeasuring personnel and comparing the original 
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and repeated measurement. If significant discrepancies are 
noted, corrective action should be taken. 

3. Replacing the integrated anthropometric measuring device 
by one using digital measurement displays, or automated 
measurement recording. 

The consequences of failing to improve the Navy's 
personnel anthropometric measuring procedure will amount to 
negating the Anthropometric Compatibility Assignment Program 
and accepting the aircraft accident rate, and ejection injury 
rate which were sufficiently undesirable to lead to the 
program's establishment. Additionally, erroneous measurement 
can adversely impact on officers' career development patterns 
and consequently on officer personnel retention. 
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APPENDIX A 



Definitions of Measurements 



1. Buttock-Heel Length -see Leg Length 

2. Buttock-Knee Length -the horizontal distance from 

the rearmost surface of the buttocks to the front of 
the kneecaps, with the subject sitting erect. 

3. Buttock-Heel Length -see Leg Length 

4. Functional Reach -the distance from the wall to 

the tip of the thumb measured with the subject's 
shoulders against the wall, his arm extended forward, 
and his index finger touching the tip of his thumb. 

5. Height -vertical distance from the 

floor to the top of the head with the subject standing 
erect in bare feet. 

6. Knee Height Sitting -the height from the footrest 

surface to the musculature just above the knee. 

7. Leg Length -the distance from the base 

of the heel to a wall against which the subject sits 
erect with his leg maximally extended forward along 
the sitting surface. 

8. Sitting Height -the height, from the sitting 

surface, to the top of the head, with the subject 
sitting erect. 

9. Stature -see height 

10. Popliteal Height Sitting -the height of the underside of 
the upper leg above the footrest surface. 
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APPENDIX B 



Anthropometric Codes as Percentiles 
of the 1964 Population 

Variable Code Measurement Equivalent 

Interval Percentile 



SITTING HEIGHT 9 

8 

7 
6 

5 

4 

3 
2 
1 
0 

FUNCTIONAL REACH 9 

8 

7 

6 

5 

4 

3 
2 
1 
0 

BUTTOCK-KNEE LENGTH 9 

8 

7 

6 

5 

4 

3 
2 
1 
0 

LEG LENGTH 9 

8 
7 

6 

5 

4 
3 
2 
1 
0 



40.0-41.0 


99+ 


39.5-39.9 


99+ 


39.0-39.4 


98-99 


38.5-38.9 


96-97 


38.0-38.4 


91-95 


35.0-37.9 


15-90 


34.0-34.9 


3-14 


33.0-33.9 


<3 


32.5-32.9 


<1 


32.0-32.4 


<1 


^27.9 


<1 


28.0-28.4 


<1 


28.5-28.9 


1-2 


29.0-29.4 


3-6 


29.5-30.4 


7-25 


30.5-30.9 


26-39 


31.0-31.4 


40-50 


31.5-32.4 


51-75 


32.5-33.9 


76-94 


^34.0 


>95 


> 28.0 


>99 


27.0-28.0 


>99 


26.5-26.9 


>98 


26.0-26.4 


97-98 


25. 5-25.9 


91-96 


25.0-25.4 


81-90 


24.0-24.9 


50-80 


23.0-23.9 


15-49 


22.0-22.9 


2-14 


^ 21.9 


<2 


49.0-50.0 


>99 


48.0-48.9 


98-99 


47.0-47.9 


95-97 


46.0-46.9 


86-94 


45.0-45.9 


71-85 


43.0-44.9 


34-70 


40.0-42.9 


3-34 


39.0-39.9 


1-2 


38.0-38.9 


<1 


36.0-37.9 


<1 
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Appendix C 



Program Listing for CAPE Model 
Used to Generate the 1964 and 1969 Data 



program GeNSAMP(FMeAS.U»-EAS.SAMPLt f INPUT t OUTPUT . 

1 Tapf1=FMEAS*TAPE2=UmEAS*TAPE J=SAPPLF t 

3 TAPrS=INPUT»TAPF6=0UTPUT) 

••• MQDiJLE A - MONTE carlo SUmULATTON ••• 

#«» MCSP READS IN Factor matrix of REOUIREO ANThpopOmetRIC 
measurements from either disk files or interactive 

TERMINAL. IT also ALLOWS USERS TO PPInT* EDIT anD SAVE 
files. the program THEN COMKUTf RANDOM MULTIVARIATE NORMAL 
«»* feature vectors, at the ENJ of The program it calls 
#»» SUPPOUTINE TRanS to perform NECESSAPY data TkANSKORmaTION 
OF anthropometric samples ANO stores The resultant 
••• measurements on disk files for cpewstation evaluation 
••• MODULE. 

COMMON/mcSm/COPP ( 78) *AnTm(2,12) *V(25) *NmEAS*KN* 

1 SOPOOT(79) .X(2S) .lU.NTMEASfiCOUNT. 

2 INAmF (4) , I TMFAC 

COMMON/MI SC/ lOI N. lOOUT. ICAPO 

NMEAS - NUMBER OF ANTHROPOMETRIC MEASUREMENTS* CURRENTLY 12 

NSAMP - NUMBER OF OPERATOR SAMPLES 

data file description: 

NAME program ref, tape# FOR/UN CONTENTS 



FMFAS 1 FOR ANT^ (2*NMEAS) *CORR (KN) 

UMEAS NTMEAS 2 UN NMEaS*KN* ANT'^*C0RP 

sample NTEMP 3 UN IN>vmE (^) *V*X 

input 5 — 

OUTPUT — 6 — 



DATA lOTN* I 00 UT*NUCREW.NTEMPfNTMEAS*NlMEAS/S* 6 t 6 * 3 * 4 . 7 / 
DATA mayFUNC/^/ 

DATA MMFAS*NTMFAS*NIMFaS/ 12 * 2 * 4 / 

1 # ICOUNT/ 0 / 

DATA AmTm/ 2 A* 0 /*CORR/ 79 * 0 / . INAME/M* 10 H / 

DATA maxOPT *MAXANTH.maxaCT/ 3 * 6 * 3 / 

♦••COMPUTE (KN) the length OF THE LINEAR ARRAY CORR(KN) 

lOOUTsO 
ISAMPsO 

•••MCSH main menu 

continue 

ICIN=5 

WRITE r<S. 290 ) 

CALL INPT ( 1 . 1 *nVALS 7 *nmfas*XVAL • lEOF) 

GO TO ( 130 , 100 * 110 . 120 . 2 M 0 . 120 ) * lEOP -1 

10 continue 

call HFLP ( 1 ) 

GO TO 100 
0 CONTINUF • 

WRIT# (^.340) 

GO TO 100 

0 continue 

IF (nmfaS.lt, 1 .OR, NMEAS. GT.MMFAS) GO TO 120 
KN=NMF AS^ (NMEAS^l ) /2 
IINPTsO 

••• GENERATE OPERATOR SAMPLE FROM mfans*STO UEV ^ CORW MATRIX 

iHO CONTINUE 
IOIN=H 

c •••anthropometric functions 

IbO CONTINUF 

WRITE <^^.300) 
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C4LL INPT (1 .1 .NVAUSZ*IOPT.XVAL. I£On 
GO TO (180, UO. 160, 170. 280* 170) . IEOF»l 
160 CONTINUP 

CALL HELP (2) 

GO TO ISO 
170 CONTINUE 

WPITP (6«3<»0> 

GO TO ISO 
180 CONTINUP 

IF ( IOPT.LT. 1 .OR. lOPT.GT.HiXANTH) GO TO 170 
C INPUT, EDIT. SAVE. PRINT, GENEkaTF OPER SAME. END 

GO TO ( 190,200,220,2^0,260.280) . lOPT 
190 CONTINUP 

c *«*inpuT anthropometric MEASUREMENTS 

CALL ACCPT (ISAMP) 

IINPT=1 
GO TO 140 
200 CONTINUE 

C ***EOIT ANTHROPOMETRIC OaTa PILE 
IF (IINPT.EO.I) GO TO 210 
WRITE (6,350) 

GO TO 140 
210 CONTINUE 
CALL EDIT 
GO TO 140 
220 CONTINUP 

C ***SAVE ANTHROPOMETRIC OATA PILE ON UmEaS 
IP (ISamP.EO.D go to 230 
WRITE (S.320) 

GO TO 140 
230 CONTINUE 
CALL SAVE 
GO TO 140 
240 CONTINUE 

C ***poiNT anthropometric OATA 
IP (IIMPT.EO.I) GO TO 250 
WRITE (6,360) 

GO TO 140 
250 CONTINUE 

CALL PRINT 
GO TO 140 
C 

C *** generate operator sample from ANTH MEAS OATA 

c 

260 continue 

IF (ISamP.ve.1) go to 270 
CALL OPSAMP (ISAMP) 

60 TO 140 
270 CONTINUP 

WRITE (6,330) 

GO TO 140 

C ***TFRMINATION OF MCSM 
280 CONTINUE 

WRITE (6,310) 

STOP 

C 

290 FORMAT ('• ENTER NUMBER OF MEASUREMENTS (1-12)? — ") 

300 FOR*^AT (•• GENERATE FUNCTIONS( 1-INPUT ;2-EOIT;3-SAVE;4-PRInTI 
- ••5-GENERaTE;6-ENO)— '•) 

310 FORMAT (•• end OF OPERATOR SAMPLE”) 

320 FORMAT (” unable TO SAVE MEASUREMENT OATA - FATAL ERRORS") 

330 FORMAT (•• unable TO GENERATE SAMPLE - FATAL ERRORS”) 

340 FORMAT (” invalid RESPONSE”) 

350 format (” ANTHROPOMETRIC DATA MUST BE- INPUT BEFORE ThE EDIT 
♦ "FUNCTION CAN BE USED”) 

E.nRM4T. (” ANTHROPOMFTRtr OATA MUST RF INPUT HFFORE ThE PRINT 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



100 



no 

c 

i?o 

no 

r 

140 

ISO 

C 

C 

C 

160 

170 

C 



♦ "FUNCTION CAN 3E USED") 

END 

SUBROUTINE INPT ( I TYPE , I VALSZ« NVALSZ » I VAL » XVAL » lEOF ) 

COMMON/MI SC/IO IN t lOOUTflCARO 

FREE FIELD INPUT ROUTINE " READS IN INPUT FROM TmE ";Eym0aR0 IN 

FREE F’lELO INPUT 

INPUT- — 

ITYPE - SPECIFIES T^E TYPE OF DATA TU HE ENTERED 
1=INTEGER 2=REAL 3=ALPHANUmERIC 
. IVALSZ - THE number OF ITEMS EXPECTED TO BE ENTERED 
OUTPUT 

NVALSZ - THE number of ITEMS ENTERED BY THE USER 
IVAL - ARRAY CONATINING INTEGER OR 

ALPHANUMERIC DATA SURPLIED BY THE USER 
XVAL - ARRAY CONTAINING REAL DATA SUPPLIED Br THE 
USER 

lEOF - END OF FILE FLAG 

1=IF BLANKS Ok A CARRIAGE CONTROL IS THE 
only INPUT 

2= IF A QUESTION Mark is entered 
3=IF INPUT IS invalid 
4=IF A $ HAS BEEN ENTERED 
DIMENSION ICHAR(80) . IVAL ( IVALSZ) .XVAL dVALSZ) 

IST0P=0 

IF (IVALSZ. GE.O) GO TO 100 
IST0P=1 

IVALSZsIABSdVALSZ) 

CONTINUE 

lOUOTEsO 

NVALSZsO 

IFOF=0 

IF (lOOUT.EO.l) GO TO 120 
DO 110 1=1. flO 
ICHAR(I)=lH 
CONTINUE 

READ IN THF input 

READ dOIN.330) ICHAR 

IF (I0IN.NE.5) WRITE (6.330) ICHAR 

IF (FOF(IOIN) .NE.l) GO TO 130 

IF0F=1 

RETURN 

continue 

I1=ICARD 

GO TO 140 

CONTINUE 

initialize the card column counters 

11 = 1 

CONTINUE 

I2=fl0- 

CONTINUE 

IF (ISTOP.EO.l. AND. IVALSZ. EO. NVALSZ) RETURN 
SFARCH FOR THE NEXT WORD 
IF (I1.LE.I2) GO TO 160 

SEARCH for THE FIRST NON-BLANK CHARACTER (START OF WORD) 

CHECK IF EXCESS HAS BEEN ENTERED 

IF (NVALSZ. LE. IVALSZ) RETURN 

TE0F=5 

RETURN 

CONTINUE 

DO 170 1=11.12 

IF (ICHaR(I) .EO.IH ) GO TO 170 

IFIRST=I 

GO TO IflO 

continue 

IF (NVALSZ. E(5.0) IE0F = 1 

CHECK IF EXCESS HAS BEEN ENTERED 
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IF (NVauSZ.GT. IVALSZ) IE0F=5 
RETURN 

IPO CONTINUF 

C SFARCH FOR LAST NON-8LANK CHARACTER (END OF wORO) 

IF ( IChaR ( IFIRST) .EO. lH?) IE0F=2 
IF (ICHAR(IFIRST).E0.1H$) IE0F=4 
IF (IF.OF.NE.O) RETURN 
IF (ITYOE.E0.4) GO TO 300 
IF {ITYPE.NE.3) go to 190 

IF <ICHAPnFlRST).NE.lH”.ANO.ICHAH(IFIRST).NE.lH') GO TO 190 
IFIRST=IFIRST*1 
IOUOTE=l 
190 CONTINUE 

00 220 I=IFIRST«I2 
IF (lOUOTE.EQ.O) GO TO 200 

IF (ICHAR(I) .EO.lH”. OP. ICHAR(I) .EO.lM') GO TO 210 
GO TO 220 

200 continue 

IF (ICHaR(I) .NE.IH .AND.ICHAR(I) .NE.lH.) GO TO 220 
210 CONTINUE 
ILAST=I-1 
GO TO 230 
220 CONTINUE 
ILAST=I? 

230 CONTINUF 

NVALSZ=MVALSZ*1 

C RESET THE FIRST CHARACTER TO RE SEARCHED FOR THE NEXT WORD 

c increment the number of words found 

•I1=ILAST*2 

IF (NVAlSZ.GT.IVALSZ) go to 150 
C CHECK IF INPUT IS alpha-numeric 

IF (ITypE.EQ.3) go to 310 
C INPUT IS NUMERIC 

C CHECK IF number is SIGNED 

ISIGN=l 

IF (ICHAR(IFIRST) .NE.lH-. AND. ICHAR(IFIRST) .rjE.lH*) 60 TO 240 
IF (ICHAR (IFIRST) .EO.lH-) ISIGN=-1 
IFIPST=IFIRST*1 
2”0 CONTINUF 

INUMsO 
IOEC=0 
IOECF=n 

DO 240 I=IFIRST»ILAST 
IF (IChaR(I) .NE. lH.) GO TO 250 
I0ECF=IDECF»1 
IDEC=1 
GO TO 260 
250 CONTINUE 

Call iconv (ichar(d .iadd.ieof) 

IF (IEOF.EQ.3) RETURN 
INUM=INUM*10*I ADD 
IF (lOEC.NE.O) IDEC=I0EC*1 
2b0 continue 

IF (IDFCE.LT. 2) 60 TO 280 
270 CONTINUF 

IEOF=3 
RETURN 

280 CONTINUE 

INUM=INUM»ISIGN 
XNUMsFLOAT (INUM) 

IF ( (ITYPE.EQ.l.AND.IDEC.LE.l) .OR. (ITYPE.EQ.2) ) GO TO 290 

IEOF=3 

RETURN 

290 CONTINUF 

IF (IDEC.NE.O) XNUM=XNUM/10.0**(IDEC-1) 

IF (ITYPF.FO.)) TVAI (NVAI S7) =TFT X (XNIIM) 
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nr>o^-^n^-o >-ooo noon 



IF ( 1 I rnt.tu.^) AVAL (NVALbif ) =XNUM 
ICAR0=ILAST*1 
GO TO 150 
300 CONTINUF 

ILAST=IVALSZ*10 
310 CONTINUF 

c INPUT 15 Alphanumeric 

DO 320 J=ltIVALSZ 
IVAL(J) slOH 
320 continue 

L=ILAST-IFIRST*1 
NVAL5Z=(L-1 ) /10*1 
IF (L.GT.10*IVALSZ) go to 270 

ENCODE (L»330.IVAL(l) ) ( ICHAR ( I ) t I =IF IPST, ILAST) 

ICARDsIl 
RETURN 
C 

330 FORMAT (80A1) 

END 

subroutine ICONV (JCHAR,II,IE0F) 

C this function converts alpha numbers to number NUMBEkS 

DIMENSION JNUM(IO) 

DATA JNUM/IHO.IHI . 1 H2 . 1 H3 t IHA , I H5 . 1H6 ♦ 1 H7» IH8« 1 H9/ 

DO 100 J=l,10 

IF (JCHaR.NE.JNUM(J) ) GO TO 100 

II=J-1 

RETURN 

100 CONTINUE 

C INVALID INPUT NOT A NUMBER 

IE0F=3 
RETURN 
END 

subroutine ACCPT (ISAMP) 

*** INPUT anthropometric DATA 
*** MEANS, STD DEV, AND CORRELATION MATRIX 



COMMON/mCSM/CORP (78) ,ANTM (2, 12) , V (25) ,NMEAS,KN, 
1 SOROOT (78) ,X (25) , IU,NTmEAS, ICOUNT, 

? INAME ( A) , ITMEAS«NIMEAS,NTEMP 

COMMON/MlSC/IOIN,IOOUT,ICARD 
DIMENSION BUFF (3) , TEMP (12) 

DATA fiUFF/3*0/ 

ISAMPsO 

ICOUNTsO 

INLINEsO 

I0IN=5 

***READ MODE - 1=INTERACTI VE 2=FILE 

00 CONTINUE 

WRITE (6*610) 

call INPT (1,1,N0INPT,IM0DE,XINPT,IE0F) 

GO TO {130,110.120.110,490,110), IEOF*l 

invalid response 

10 WRITE (6,650) 

GO TO 100 

HELP 

20 CALL HELP (4) 

GO TO 100 

30 IF (IMODE.EQ,!) GO TO 200 

IF (IM0DE,NE,2) GO TO 110 



*** file Type: i=formatted 2=unformatteo • 



140 



CONTINUF 



67 



WRITE (6.620) 

CALL INPT (I.l.NJOlNPT.IFFORM.xINPT.IEOn 
GO TO (170,150.160,150,490,150). IE0F*I 
150 WRITE (6,650) 

GO TO 140 

160 CALL HELP (5) 

GO TO 140 

170 IF (IFFORm.EO.I) go TO 190 
IF (IFFORH.NE.2) go to 150 
C 

c •** unformatted read 

REWIND NTMcas 

c ***rpad The number of anthropometric measurements, 
c ***correlation components 

read (ntmeas) nmeas.kn, iname 

IF (EOF (NTmEaS) ,NE. I ) GO TO 180 
WRITE (6.670) 

GO TO 490 

C ***RFAD ThF mean and STANDARD DEVIATION VALUES 
180 RFAD (NTMEAS) ( ( ANTm ( 10 , JO ) , 1 0= 1 »2 ) » JO=I ,NME AS) 

c ***reaD The correlation matrix values 

READ (NTMEAS) (COPR ( 10) , 10=1 ,KN) 

REWIND NTMEAS 
WRITE (6,690) 

I5AMP=l 

RETURN 

C *** SET INPUT LOGICAL UNIT TO I FOR FILE FmEAS 
190 IOIN=I 
200 CONTINUE 
C 

C ***ACCEPT operator SAMPLE FILE NAME 
C 

call NAMEFL (IERROR) 

ICOUMTalCOUNT*IERROR 

INLINE=INLINE*1 

c *** inputs means and standard deviations into array ANTm 
IF (IMOOE.EO.l) write (6,500) 

IF (IMOnE.NE.l) WRITE (6,510) INAME 
MAXINPT=2 

DO 320 IROw=I,NMEAS 
210 CONTINUE 

IF UMODE.EO.l) WRITE (6,530) IROW 
220 NAMT=0 
230 CONTINUE 

ITOTINsmaxinPT-NAMT 

IC0L=NAMT*I 

CALL INPT (2,IT0TIN,N0INPT,IVAL,ANTM(IC0L,IR0W) ,ItOF) 

IF (IEOF.EQ.5.0P.IEOF.EO.O) NAMT=NamT*NOINPT 

IF (IMOOE.EO.l) GO TO 240 

INLINE=INLINE*I 

IF (NAmT.EO.MAXINPT) go TO 290 
IF (lEOF.EO.l) GO TO 480 
WRITE (6.540) IROW 
IF (IEOF.EO.S) go TO 270 
WRITE (6,700) INLINE 
IC0UNT=IC0UNT*1 
60 TO 320 

C*** INTERACTIVE lEOF PROCESSING 
240 CONTINUE 

GOTO (280,210,250,260,490,270), IEOF*I 
250 CONTINUE 

CALL HELP (6) 

GO TO 210 
260 CONTINUE 

WRITE (6»650) 

GO. TO 210 
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270 CONTINUE 

IOLISTsNOInPT-ITOTIN 
W°ITE (6«600) lOLIST 
IF (IMODE.NE.l) ICOUNT=ICOUNT*l 
GO TO 2<?0 
280 CONTINUE 

IF . (NAMT.GE.2) go to 2<?0 
WRITE (fttSSO) 

GO TO 230 
290 CONTINUE 

IF (NAMT.GT.2) NAMT=2 
IRERP=0 

DO 300 ICOl.= l«NAMT 

CALL PFASON ( 1 .ICOLt IPOW»ANTM( 1 .IROW) .ANTM(2«IR0 W) . IFaTAL) 
IRERR=IPERP*IFATAL 
300 CONTINUF 

IF (IREPP.EO.O) GO TO 320 
IF (IMOOE.NE.l) GO TO 310 
WRITE (6.630) IROW 
GO TO 220 

310 ICOUNT=ICOUNT*IRERR 
320 CONTINUE 
C 

C ANTHROPOMETRIC CORRELATION MATRIX VALUES ARE ACCEPTED 

C *** AND PROCESSED. 

C 

IF (IMOOE.EO.l) WRITE (6.560) 

IF (IMODE.NE.l) WPITE (6.520) 

00 A60 lROw=l.NMEAS 
330 CONTINUF 
NAMT=0 

MAXINPT=NMEAS*1-IR0W 
3A0 CONTINUE 

IF ( ImoDE.EQ. 1 .AND. IROW. EO.NMEAS) WRITE (6.580) 1 80>’ . MAX INPT .NME AS 
IF (IMOOE.EO.l. AND. IROW. NE.NMEAS) WRITE (6.570) I wOw .MAX INPT . IPO- . 
INMEAS 

350 CONTINUF 

itotin=m4xinpt-namt 

ISTART=IROW*NAMT 

CALL INPT (2.IT0TIN.N0INPT.IVAL.TEMP(ISTART) ,IEOF) 

IF ( lEOF.EO.O.OR. IEOF.EO.5) NAMT=NAMT»N0 INPT 

IF (IMonE.EQ.l) GO TO 360 

INLINF=INLINE*1 

IF (MAmT.EO.MAXINPT) go to 410 
IF (lEOF.EO.O) GO TO 350 

IF (lEOF.EO.l) WRITE (6.710) 

IF (lEOF.EQ.l) GO TO 470 

IF (IEOF.EO.5) GO TO 390 

ICOUNT=ICOUNT*l 
WPITE (6.700) INLINE 
GO TO 460 

C *** interactive input error PROCESSING 
360 CONTINUF 

GO TO (400.380.370.380.490.390). IEOF»l 
370 CONTINUE 

CALL HELP (7) 

GO TO 340 
380 CONTINUE 

WRITE (6.650) 

GO TO 340 
390 CONTINUE 

IOLISTsNOINPT-ITOTIN 

WRITE (6.600) lOLIST 

IF (IMODE.NE.l) ICOUNT=ICOUNT*l 

NAMT=MAXINPT 

GO TO 410 
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400 

410 

420 

430 

440 

450 

4A0 

470 

480 

490 

c 

500 

510 

520 

530 

540 

550 

560 

570 

5>i0 

590 

600 

610 

6?0 

630 

640 

6S0 

660 

670 

640 

690 

700 

710 



CONTINUE 

IF (NOTNPT.EO.ITOTIN) 60 TO 410 ’ 

IT0TIN=m4XINPT-NAMT 
WPITE (6*590) ITOTIN.IPOW 
60 TO 350 
CONTINUE 

ILAST=I«0W*N4MT-1 

irerr=o 

DO 420 IC0l=IB0W.ILAST 

CALL REASON (2*IC0L*IB0W*TEMP(IC0L> *0*IE4TAL) 

IRERR=IRERR*IEATAL 

CONTINUE 

IE (IRERR.EO.O) 60 TO 440 
IE (IMOnE.EO.2) 60 TO 430 
WRITE (6*640) IROW 
GO TO 3"<0 

ICOUNT=ICOUNT*IRERR 
60 TO 460 
CONTINUE 

DO 450 TCOL=IROW*NmEAS 
N=NDX ( IROW, ICOL) 

CORR(N)=TEmP(TCOL) 

CONTINUE 

CONTINUE 

CONTINUE 

IE (ICOUNT.EQ.O) ISAMP=1 
IE (IMODE.NE.l) WRITE (6*680) ICOUNT 
60 TO 490 
WPITE (6*660) 

return 

CONTINUE 

RETURN 

FORMAT (•• ENTER PAIRS OF MEANS AND STANDARD DEVIATIONS — ") 
format (//" DIAGNOSTICS OF FORMATTED READ ”/•• OPER OESC: "*4A10/) 
FORMAT (//•• CORRELATION/ANTHROROmETRIC MEASUREMENTS!") 

FORMAT (" measurement - "*I2*" ? ") 

FORMAT (" measurement -"*12) 

FORMAT (" ENTER standard DEVIATION — ") 

FORMAT (" ENTER ANTHROPOMETRIC CORRELATION MATRIX DATA IN A ROWWIS 
IE UPPER "/" TRIAN6ULAR FORM."/" NOTE! ALL DIA60NAL ELEMENTS MUST"* 

♦ " FOUAL 1.000") 

FORMAT (5h row *I2*"("*I2*" ELEMENTS* COLUMNS " * 12* "-" * 12 *") " ) 
FORMAT (5H POW *12. "(",12*" ELEMENT* COLUMN "*I2*")") 

FORMAT (2lH ENTER THE REMAININO *I2*1HH ELEMENTS EOR ROw *12) 
FORMAT (" The last "*I?*" INPUTS HAVE BEEN I6NORED") 

FORMAT (" input MOOE(l=INTERACTIVt;?=FlLE)— ") 

FORMAT (" READ MODE ( 1=F0RMATTED!2=UNF0RmaTTED) — "> 
format (" re-enter data for "*I2) 

FORMAT (" RE-ENTER DATA FOR ROW "*I3) 

FORMAT (" invalid RESPONSE") 

format (" UNEXPECTED EOF ON FMEaS FOR MEAS "*I2) 

FORMAT (" FILE UMEAS IS NOT ATTACHED TO ThE PR06RAM"/ 

♦ " unformatted read ON THIS FILE IS ABORTED "*//) 

FORMAT (" TOTAL ERRORS ON INPUT JATA ",I5) 

FORMAT (" MEAN*STD DFv 6 CORRELATION DATA READ") 

FORMAT (" invalid DATA AT LINE "*I5) 
format (" unexpected end OF FILE") 

END 

SUfaPOUTINE EDIT 

COMMON/MCSM/CORR (78) *ANTM(2*12) *V(25) *NM£AS*KN* 

1 SOROOT (78) *X (25) ♦ IU.NTmEAS* ICOUNT * 

2 INamE(4) ,ITMEAS*NImEAS*NTEMP 
COMMON/MISC/IOIN, lOOUT.ICARO 
DIMENSION BUFF(3) 

DATA BUFF/ 3*0/ 
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c 

C “♦* MOOIflCATIONS ON MEflNStSTANOAPO DEVIATIONS* 

C *** AND ANTHROPOMETRIC CORRELATIONS 

C 

IOIN=5 

100 CONTINUE 

WRITE (6.350) 

CALL INPT (l.l.NVALSZ.IOPT.XVAL.IEOE) 

GO TO (130,100.110.120,340.120). IEOF*l 
no CONTINUE 

CALL help (8) 

GO TO 100 
120 CONTINUE 

WRITE (6,360) 

GO TO 100 
130 CONTINUE 

IE (IOPT.LT.l.OR.IOPT.GT.5) GO TO 120 
60 TO (140,140.240,330,340). lOPT 
C 

C *** ICOL=l EOR MEAN PROCESSING 

C *** IC0L=2 EC^ standard deviation PROCESSING 

C 

IwO CONTINUE 
IC0L=I0PT 
150 CONTINUE 
M4 XInPt=2 
NAMTsO 

160 CONTINUE 

IE (ICOL.EO.1) write (6.380) 

IE (ICOL.EO.2) write (6.3R0) 

170 CONTINUE 

I VALSZ=maXINPT-NAMT 
J=NAMT+1 

call INPT (P.IVAL'sZ.NVALSZ.IVAL.SUEElJ) .lEOE) 

GO TO (220,180.200.190,100,210) , IEOE*l 
160 CONTINUE 

WRITE (6.370) 

GO TO 160 
190 CONTINUE 

WRITE (6.360) 

GO TO 160 
200 CONTINUE 

IE (ICOL.EO.1) call help (11) 

IE (ICOL.EO.2) call help (12) 

GO TO 160 
210 CONTINUE 

IOLISTsNVALSZ-IVALSZ 
WRITE (6.400) lOLlST 
220 CONTINUE 

IPOW=RilEE(l) 

NAMT=NVALSZ*NAMT 
IE (NAMT.GE.MAXINPT) go to 230 
IE (ICnL.EO.l) WRITE (6,420) IROW 
IE (ICOL.EQ.2) write (6.430) IROW 
GO TO 170 
230 CONTINUE 

CALL REASON (1.ICOL,IROW,8UEE(2) .O.IEATAL) 

IE (IEaTAL.EO.1) go TO 150 
ANTM(ICOL.IROW)s5UEE(2) 

GO TO 100 
C 

C *** processing EOR ANTHROPOMETRIC CORRELATIONS 
C 

240 continue 
MAXINPT=3 
NAMT=0 
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aso continue 

WPITE (*>.410) 

260 CONTINUE 

IVaLSZ=MAXlNPT*N4MT 

j=NAmT*1 

CALL INPT (2.IVALS7.NVALSZ»IVAL«3UFF(J) .lEOF) 
GO TO (310,270.290,280,100.300) ♦ IEOF*l 
270 CONTINUE 

WPITE (6,370) 

GO TO 250 

280 continue 

WRITE (6.360) 

GO TO 250 
290 continue 

CALL HELP (13) 

GO TO 250 
300 CONTINUE 

I0LIST=NVALSZ-IVAL5Z 
WRITE (6.400) lOLIST 
GO TO 320 
310 continue 

ICOL=0UFE(1) 

IPOW=BUFF(2) 

NAMTsNAMT.NVALSZ 
IF (NAMT.GE.MAXINPT) go to 320 
IF (NAmt.EO.1) write (6.440) ICOL 
IF (NAMT.EQ.2) write (6.450) ICOL.IROW 
GO TO ?60 
320 CONTINUE 

CALL reason (3. ICOL. IR0W.8UFF(3) .O.IFATAL) 

IF (IFaTAL.EQ.I ) GO TO 240 

CALL RFASOn (2. ICOL. IP0W.8UFF(3) .O.IFATAL) 

IF (IFATAL.EQ.I) GO TO 240 

C ***LOAO The value of buff (3) INTO T-<E APPROPRIATE 
c ***P0SITI0N within the CORR(N) array 





N=NOX ( IROW 


tICOL) 




CORR (N) 


=RUFF(3) 




GO TO 100 




c 

c 


**« print 


OPTION 


L 

330 


CONTINUE 






CALL PRINT 






GO TO 100 




C 


«»* END 


OF 


EDIT mode 


340 


CONTINUE 




r- 


retupn 






L 

350 


FORMAT 


(•* 


EDIT MODE (1 -MEANS ;2-STU:3-ANTh. CORR? 4-PRInT I 5-ENO) 


340 


FORMAT 


(•* 


invalid RESPONSE") 


370 


format 


(•* 


NO INFORMATION ENCOUNTERED") 


3H0 


format 


(♦* 


INPUT INDEX AND MEAN VALUE") 


390 


format 


(•• 


INPUT INDEX AND STANDARD DEVIATION") 


400 


format 


(*• 


THE LAST". 12." INPUTS HAVE BEEN IGNORED") 


410 


FORMAT 


(** 


INPUT COLUMN . ROW AND CORRELATION") 


4^0 


FORMAT 


(** 


INPUT MEAN VALUE FOR INDEX ".I3) 


430 


format 


(** 


INPUT standard deviation for index ".13) 


440 


format 


(•* 


INPUT ROW AND CORRELATION FOR COLUMN ".12) 


450 


format 


(•* 


INPUT CORRELATION FOR ROW ".12." COLUMN ".12) 




END 








SUBROUTINE 


NAMEFL (IFPROR) 




COMMON/MCSH/CORR (70) .ANTM(2.12) .V(25) .NMEAS.KN. 



1 SOROOT(7fl) .X(25 ) .IU.NTmEAS.ICOUNT. 

2 INamE (4) . ITMEAS.NIMEAS.NTFMP 
COMMON/mISC/IOIN.IOOUT.ICAPO 
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C ***PROMPTS USER FOR NAMf OF SAMPLE FILE 
C 

DIMENSION IRUFF (4) ,XVAL(4) 

DATA IRUFF/4*10H / 

IFRRORsO 

C 

C ***INITIalIZE CONTENTS OF I9UFF(4) JITH 
c ***PLANK characters 
c 

DO 100 

IRUFF( J)=10H 
100 CONTINUE 

c 

no CONTINUE 

IF (lOIN.EO.S) WRITE (6.180) 

CALL INPT (4,4.NVALSZ.I8UFF.XVAL»IE0F) 

IF (IE0F.LT.2.0R.IE0F.E0.5) 00 TO ISO 
IF (I0IN.NE.5) GO TO 140 
60 TO (120,130.170), lEOF-l 
120 CONTINUE 

CALL help (14) 

GO TO no 
130 CONTINUE 

WRITE (6.190) 

GO TO no 
140 CONTINUE 

WRITE (6.190) 

IFRRORsl 

C ***LOAO contents OF I8UFF INTO INAM£ ARRAY 
150 CONTINUE 

DO 160 J=l,4 
INAMF(j)=I3UFF(J) 

160 CONTINUE 
170 RETURN 
C 

1«0 FORMAT (■• OPERATOR SAMPLE OESCRIPTION(40 CHA)< MAX) '•) 

190 FORMAT ('• invalid RESPONSE") 

END 

SUBROUTINE OPSAMP (ISAmP) 

c ***generate operator sample 

COMmON/mCSM/CORP (78) ,ANTM(2.12) ,V(25) .NMEAS.KN. 

1 SOROOT (78) .X (25) . lU.NTME AS. ICOUNT. 

2 INAME (4) , I TMEAS.NIMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARO 
DIMENSION IRAN(9) 

DATA ( IRAN( I ) .1 = 1 .9) /33 333333 .55555555 . 7777777 7 . 22222222 1 . 

1 444444443.666666669.886888689,999999999.123456789/ 

C 

C ««*LINEAR array C0RR(78) IS LOADED INTO LINEAR 
C »«*DUMMY array SOROOT(78) BEFORE MATRIX TRANSFOR- 
C »»*MATI0NS. all MATRIX TRANSFORMATIONS ARE 
C ««*ppreoRmED on SOROOT (78) TO PRESERVE C0RR(78) 

C ***AS implemented within the editing PROCESS. 

c 

I0IN=5 

DO 100 K=1.KN 
S0R00T(K)=C0RR(K) 

100 CONTINUE 

C perform MATRIX TRANSEORMATIONS 

call SOPMX (NERROR) 

IF (NF.RROR.EO.1 ) return 
c *** accept operator sample size 
no continue 

WRITE (6.320) 

CALL INPT (l.l.NVALSZ.NSAMP.XVAL.IEOF) 

GO TO (150. 1?0. 130. 140. PRO. 140) . IFOF.l ' 
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120 CONTINUfi: 

WPITE (6.330) 

GO TO 110 
130 continue: 

CALL HELP (IS) 

GO TO 110 
140 CONTINUE 

WPITE (6.340) 

GO TO 110 
150 CONTINUE 

IE (NSAmP.GT.O) GO TO 160 
WPITE (6.350) 

GO TO 110 
C 

C •**ASK USEP EOP machine TYPE (32 OP 60 BIT) 

C 

160 CONTINUE 

WPITP (6.300) 

CALL INPT (l.l.NVALSZ.MTYPE.XVAL.IEOE) 

GOTO (190.160.170.180.290.180). IEOF»l 
170 CONTINUE 

CALL HELP (20) 

GO TO 160 
180 CONTINUE 

WPITE (6.340) 

GO TO 160 
190 CONTINUE 

IE («TYPE.NE.32. ANO.MTYPE.NE.60) GO TO 180 

200 continue 

WPITE (6.310) 

C *** OPTION EOP accepting A PaNOOm SEED. 

CALL TNPT ( 1 . I .NVALSZ. IS.XVAL. lEOF) 

GO TO (240.210.220.230.290.230). IEOF*l 
210 CONTINUE 

WPITE (6.330) 

GO TO 200 
220 CONTINUE 

CALL HFLP (21) 

GO TO 200 
230 CONTINUE 

WPITE (6.340) 

GO TO 200 
240 CONTINUE 

IF (IS.EO.O) GO TO 250 

IE (IS.LT.I.OP.IS.GT.9) go to 230 

IU=IRAN(IS) 

GO TO 260 
250 CONTINUE 

IU=671I1I33 
260 CONTINUE 

REWIND NTEMP 

C ***WPITE 40 CHARACTER NAME FIELD INTO 
C «*»FILE sample (NTEMP) 

WPITE (NTEmP) INAME.NMEA5 
C 

DO 280 I=1.NSAMP 

C *•• SIMULATE A PSEUDO SAMPLE CASE. 

CALL SIMSU3 (MTYPE) 

DO 270 K=1.NM£AS 
V(K)=V(K)*ANTM(2.K) .ANTMd.K) 

270 CONTINUE 

WPITE (NTEmP) (V( 10) .I0=1.NMEAS) 

280 CONTINUE 

C «**SET ISamP to 1 TO INDICATE OPERATOR SAMPLE HAS BEtN GtNEPATEO 

ISAMPri 

END Fit E.. NTEMP 
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REWIND 


ntemp 










REWIND 


ntmeas 










WRITE 


(6.360) 








?90 


CONTINUE 












RETURN 












300 


format 


(*• 


ENTER 


machine word size (32 


OR 60 bIT) - 


-.1) 


310 


FORMAT 


(*• 


ENTER 


RANDOM NO SEED (0 IF 


NONE OH 1-9)- 


• II) 


3?0 


FORMAT 


(** 


INPUT 


SAMPLE SIZE — ") 






330 


FORMAT 


(** 


NO information ENCOUNTERED" 


) 




3^0 


FORMAT 


(*• 


invalid RESPONSE") 






350 


FORMAT 


(*• 


sample 


SIZE MUST BE greater 


THAN 0") 




3^0 


FORMAT 


(•' 


SAMPLE 


DATA SAVED ON FILE SAMPLE OR TaPL 


3.*») 




END 














SUBROUTINE 


PRINT 









c 

c ***OUTPUT OF SIMULATED OPERATORS 
C 

C **-*NOTE: ALL OPERATOR MEASUREMENTS 

C ***ARE STORED ON FILE SAMPLE. PRESENTLY EACH RECORD 

c *** contains nmea'' anthropometric measurements 
c 

C0MM0N/MCSM/C0RR(7fl) ♦ANTM(2»12) *Y(25) .NMEAS.KN. 

1 S0R00T(78) »X(25) . lU.NTME AS » I COUNT , 

2 INA^E (A) .ITMEAS.NImEAS.NTEMP 
COMMON/MISC/IOIN. lOOUT.ICARD 
INTEGER FMT 

DIMENSION FMT(f>) .NUM(1?) .TEMP (12) .INDEX (2) 

DATA FMT/AHC .AH .AH(6X).AH .AH (F6. .AH3) ) / 

DATA NU«/1H1 .IH2.1H3.1HA.1H5.1H6.1H7.1HP.1HS'.2H10.2H11.2H12/ 
ISTOREsO 
IOFLAG=0 
C 

C *** MENU CHOICE OF OUTPUTS 
C 

100 CONTINUE 

WRITE (6.630) 

CALL INPT (l.l.NVALSZ.IOPT.XVAL.IEOF) 

60 TO (120.110.110.110.390.110). IE0F*1 
110 CONTINUE 

WRITE (6.650) 

GO TO 100 
120 CONTINUE 

IF ( lOPT.LT. 1 .OR. lOPT.GT.A) GO TO 110 
IF (lOPT.EO.l) GO TO 130 

IF (I0PT.E3.2) GO TO 150 

IF (lOPT.EO.A) GO TO 390 

GO TO 190 
C 

C ♦** printout of MEANS/STANDARD DEVIATIONS 

C 

130 CONTINUE 

WRITE (6.580) 

DO lAO iROWsl.NMEAS 

WRITE (6.6A0) IROW. (ANTm(JO.IROW) .J 0=1.2) 
lAO CONTINUE 

WRITE (6.620) 

GO TO 100 
C 

C *** RRINTOUT OF ANTHROPOMETRIC CORRELATIONS 

C 

150 CONTINUE 

WRITE (6.590) 

WRITE (6.660) (NUM( 10) . 10=1 .NMEAS) 

WRITE (6.600) 

00 180 IROWsl.NMFAS 
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00 160 TC0L=l«NMEftS 
N=N0X ( IPOW, ICOL) 

TEMP(IC0L)=C0RR(N) 

160 CONTINUE 

IF (IROW.GT.l) GO TO 170 

WRITE (6.610) (TEMP ( ICOL) . ICOL = I-^Ow.NMEAS) 

GO TO 1«0 

170 FMT(2)=NUM( IROW-1) 

FMT(4) =NUM(NMEAS*1-IR0W) 

WRITE (6.FHT) (TEMP(ICOL) .IC0L=IR0W.NMEAS) 

160 CONTINUE 

WRITE (6.620) 

GO TO 100 
C 

c *** PRINT OF anthropometric MEASUREMENTS 

c 

190 CONTINUE 

WPITF (6.430) 

IF (lOFLAG.LT.l) WRITE (6.440) 
lOFLAGsl 
C 

C **«impuT CARRIAGE" RETURN — ENTIRE SAMPLE LIST 
c ***InPUT single number — single OPERATOR 

C **«inPUT TWO numbers — SEGMENT OF OPERATOR LIST 

C 

CALL inrt (1.2. nvalsz. index. XVAL.IEOF) 

GO TO (230.270.200.210.100.220) . IEOF*l 
200 CONTINUE 

CALL HELP (19) 
lOFLAGsO 
GO TO 190 
210 CONTINUE 

WRITf' (6.420) 

GO TO 190 
220 CONTINUE 

I0LIST=NVAlSZ-2 
WRITE (6.450) IOLIST 
230 continue 
c **«imoex validation 

IF (INOFX(l) .GT.O) GO TO 240 
WRITE (6.480) 

GO TO 190 
2U0 CONTINUE 
C 

c ***TO DETERMINE WHETHER FILES WILL 9£ REWOUND 
C 

IF (ISTORE. GT.O. AND. INOEX(l) .GT.ISTORE) GO TO 250 
NRSKIP=INDEX(1)-1 
GO TO 260 
250 CONTINUE 

NPSKIP = INDEX(1) -I STORE-1 
260 CONTINUE 

ISTAPT=IN0EX(1) 

C ***IF TWO INPUTS ENCOUNTERED - OUTPUT LIMITS 
C ***ip ONE INPUTS ENCOUNTERED - OUTPUT SINGLE OPERATOR 
IF (NVALSZ.EO.1) go to 280 
IF (NVALSZ.GE.2) go to 290 
270 CONTINUE 
ISTOREsO 
ISTARTsl 
IEND=999 
NRSKIPsO 
GO TO 320 
200 CONTINUE 

IEND=ISTART 
GO TO S10 
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2<i0 CONTINUE 

IF IlNOFXd) .LE.INDEX(2) ) GO TO 300 
WPITE (6«490) 

GO TO 190 
300 CONTINUE 

IEN0=IN0EX(2) 

310 CONTINUE 

IF (INOEXm .GT.ISTORE.flNO.ISTORE.NE.O) GO TO 330 
320 CONTINUE 
C ***REWINO FILE • 

REWIND NTEMP 

READ (NTEMP) INAMEtNMEAS 
IF (EOF(NTEMP) .EO.l) GO TO 410 
330 CONTINUE 

IR (NPSKIP.LT.I) GO TO 350 

c ***wiNo File forward 

DO 340 1=1 ,NRSKIP 

READ (NTEMP) (V(K) »K=1 ,NMEAS) 

IF (EOF (NTE»^P) .EO.l ) GO TO 400 
340 CONTINUE 
3S0 CONTINUR 
C ***OUTPUT header 

WRITE (6«530) (INAME(IO) .10=1.4) 

IF ( ISTaRT.EO. lEND) WRITE (6.540) lENO 

IF (ISTAPT.NE.IENO. AND. IEN0.NE.9V9) wPiTE (6.550) ISTART.IENO 
WPITF (6.500) 

WRITE (6.660) (NUM(IO) .IO=l.NMEAS) 

WRITE (6.600) 

C 

C ***OUTPUT DATA FROM FILE 
C 

ISKIP=0 

NFOUNO=0 

00 360 I=ISTART. lENO 
C ***SKIP LINE REATURE 

IF (ISKIP.RO.l) WRITE (6.460) 

ISKIP=0 

IF (ISTART.NE.IENO.ANn.I.NE.IENO.ANO.(I/5)*5.EO.I) ISMP=1 
READ (NTEMO) (V(K) .K=l . nMEAS) 

IF (EOF(NTRMP) .EO.l) GO TO 370 

NFOUNO=NFOUND*l 

ISTORE=I 

WPITE (6.510) (V(K) .K=l .NMEAS) 

360 CONTINUE 
GO TO 3fl0 
370 CONTINUE 

IF (IENn.NE.999.AN0.ISTART.NE.IEN0) WRITE (6.470) lENO 
ISTOPE=I 

IF (NFOUNO.lt, 1 ) WRITE (6.560) 

380 CONTINUE 

WPITE (6.520) 

GO TO 190 
390 continue 

rewind ntemp 

RFTURN 

400 CONTINUE 

I0lIST=I-2*IST0RE 
WPITE (6.560) 

GO TO 190 
410 CONTINUE 

WPITE (6.570) 

GO TO 100 
C 

420 format (" invalid RESPONSE") 

430 format (» OUTPUT MODE: — ") 

uuO FORMAT (" FOR AIL OPFRATOR DATA (ENTER CARR I AGE RETURN) »./ , 
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^50 

460 

470 

4fl0 

490 

500 

510 

520 

530 

540 

550 

560 

570 

5M0 

590 

600 

610 

6?0 

630 

640 

650 

660 



C 

C 

C 

C 

C 



100 



no 



c 

c 

c 



c 

c 

r. 



1" FOR INDIVIDUAL operator DATA (ENTER OPERATOR NUHBER)”*/ 

2" FOR A SEQUENCE OF OPERATORS (ENTER RANOE OF OPERATOR NUMBERS)"/ 

* " END (ENTER $)") 

FORMAT (•• THE LAST"»I4," INPUTS HAVE BEEN IGNORED") 
format (1H ) 

format (•• error in upper index — ".14) 

FORMAT (" INDEX number CANNOT BE LESS THAN OH EQUAL TO 0") 

FORMAT (" FIRST INDEX MUST BE LESS THAN SECOND INUEX") 

FORMAT (27x, "OPERATOR SAMPLE" ./. 71 ( 1 h-) . ) 

FORMAT (11F6.2.F5.2) 

FORMAT (71 (IH-)) 

FORMAT (1X.4A10.10X.A10) 

format ( lOX. "OPERATOR NO". 14) 

format ( lOx. "OPERATOR NOS". 14." — ".14) 

format (" INDEX NUMBER EXCEEDS RECORDS IN OPERATOR FILE") 
format (" FILE SAMPLE IS NOT ATTACHED TO THE PROGRAM"/ 

♦ " operator OUTPUT REQUEST ABORTED ".//) 

FORMAT (/.21X. "MEANS AND STANDARD DEV I AT IONS" ./ . 1 X . 70 ( IH- ) ./ .2X . "M 



1EASUREMENT".6X. "MEANS". 6X. "STD DEV"./.lX.71 (1H-) ) 

FORMAT (/.7X. "CORRELATION MATRIX FOR ANTHOPOMETR IC MEASUREMENTS"./ 
1.1X.70(1H-) ) 

FORMAT (1X.70(1H-)) 

format (12F6.3) 
format (1X.70(1H-)) 

roRMAT (" PRINT MOOE( 1-MEAN/STD;2-C0RR MATRXI3-0PER mEAS ;4-END) -" ) 
format (I9.5X.2F12.2)) 

FORMAT (" INVALID RESPONSE") 

FORMAT (2X.11 (A2.4X) .A2) 

END 

subroutine RAND2 (RNORM.MTYPE.OUAN. I J) 



•••generates a normal DISTRIBUTION -kITH MEAN OF ZERO AND 
•••STANOAPO DEVIATION EQUAL TO ONE.... 

•••METHODOLOGY MARSAGLIA - BRAY 



C0MM0N/mCSm/C0RR(7R) .amTM(2.12) .v(2S) .NmEAS.KN. 

1 S0R00T(78) .X(25) . lU.NTME AS . ICOUNT . 

2 INAMEI4) .ITMEAS.NIMEAS.NTEMP 
DIMENSION 0UAN(2) 

COMMON/MISC/IOIN.IOOUT.ICARD 

IJ=IJ*1 

IF (IJ.FQ.2) go to no 

CONTINUE 

X 1=2.0 •UNFORM (MTYPE) -1 .0 
Yl=2.0^UNF0RM(MTYPE) -1.0 
S=X1*X1*Y1^Y1 
IF (S.GE.1.0) 60 TO 100 
S=SORT(-2.0^ALOG(S)/S) 

0UAN(1)=X1*S 
OUAN(2) sYl^S 
CONTINUE 
RN0RM=0UAN(IJ) 

IF (IJ.EQ.2) IJ=0 

RETURN 

END 

SUBROUTINE REASON ( ICALL . ICOL ♦ IROW . TSTORE 1 . TST0RE2 . IFAT^L ) 

••• CHECKS FOR reasonableness OF INPUT VALUES 

COMMON/MCSM/COPR (7B) .ANTM(2.12) .V(25) .NmEAS.KN. 

1 SOROOT (7B) .X (25) . lU.NTwf AS. ICOUNT. 

? INAmE (4) .ITMEAS.NIMEAS.NTEMP 

COMMON/mISC/IOIN.IOOUT.ICARD 



ICALL = 1 MEANS AND STANDARD DEVIATIONS 
TCALI = ? CORREI ATTON MATRT.X. VFPTFTCATTON 
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C ICALL = 3 INDEX VEi^IFICATION 
C 

IFAT4L=0 

lAFLAGsO 

IF (ICALL.lt. 0) IAFLAG=1 
ICALL=IABS(ICALL) 

IF (Icall.gt. 1 ) go to no 

C 

C •** mfan values must be within a three standard deviation 

C *** neighborhood of The ACCEPTEO mean values 

c 

IWRITE=1 

IF (IROW.lt. 1. OR. IROW.GT.NMEAS) GO TO 120 
IF (ICOL.EO.2) go to 100 

return 

c 

c *** STANDARD deviations MAY NOT BE NEGATIVE 

c *** AND The means must be greater Than 

c *** The standard deviation. 

c 

100 CONTINUE 
IWRITE=2 

IF (TSTflRE2.LT.O.) GO TO 120 
IWRITE=3 

IF (TST0RE1.LE.TST0RE2) GO TO 120 
RETURN 
C 

c *** CHECKS validity of index numbers 

C 

110 CONTINUE 
IWRITE=1 

IF (IROW.LE.O. OR. IROW.GT.NMEAS) GO TO 130 
IWRITE=? 

IF (ICOL.LE.O.OP.ICOL.GT.NMEAS) go to 130 

c ***is The number within the upper triangular region? 

IWRITE=3 

IF (IPOW.GT.ICOL) GO TO 130 
IF (ICalL.EO.3) return 
C 

c *** all diagonal elements must EUUAL 1.000 

C *** matrix elements May NOT BE LESS THAN (-1) OR GREATER 

C •*» Than one. 

C 

IWRITEsA 

IF (ICOL.EQ.IROW. AND. TSTOREl.NE. 1.000) GO TO 130 
IWRITFrS 

IF (TSTORE1.lt. -1.000) GO TO 130 
IWRITE=6 

IF (TSTOREl.GT. 1.000) GO TO 130 

return 

C 

c ***branch to appropriate error message 

c 

120 CONTINUF 

IF (IWRITE.NE.l .AND.IOIN.NE.5) WRITE (6.280) IROW 
GO TO (190,140,150), IWRITE 
130 CONTINUE 

GOTO (190,200,210,160,170,180), IWRITE 
140 CONTINUE 

WRITE (6,230) 

GO TO 220 
150 CONTINUE 

WRITE (6,240) 

GO TO 220 
160 CONTINUE 

WPITF (6.?S0) .T.ROW.irOl 



79 





60 TO 


220 






WO 


CONTINUE 








WPITF 


(6.260) 


IROW.ICOL 






GO TO 


220 






IdO 


CONTINUE 








WRITE 


(6.270) 


IROW, ICOL 






GO TO 


220 




• 


190 


CONTINUE 








WRITE 


(6.290) 


IROW.NMEAS 






GO TO 


220 






200 


CONTINUE 








WPITE 


(6.300) 


ICOL.NMEAS 






GO TO 


220 






210 


CONTINUE 








WPITE 


(6.310) 


IROW.ICOL 






GO TO 


220 






220 


CONTINUE 








IFATAL 


:^1 








return 








V 

230 


FORMAT 


(" standard deviations less Than ZERO") 


240 


format 


(" STfc.JDARD DEVIATION GREATER 


Than mean value.") 


250 


FORMAT 


(" DIAGONAL ELEMENT IN ROw ". 


12." COLUMN ".12. 




♦ "IS 


NOT EQUAL TO 1.0") 




260 


format 


(" ENTRY IN COLUMN ".12." ROW 


".12. 




♦ "IS 


LESS than -1.0") 




270 


format 


(" entry in COLUMN ".12." ROW 


".12. 




♦ " IS 


greater 


than .1.0") 




280 


format 


(IIO) 






290 


format 


(" ROW 


".12." IS LESS than zero OR GREATER THAN ". 


300 


FORMAT 


(" COLUMN ".12." IS LESS THAN 


ZEPO OR GREATER THAN 


310 


FORMAT 


(*• ROW 


".12." COLUMN ".12." IS OUTSIDE OF UPPER". 



♦ •' triangular region of correlation Matrix.") 

END 

SUBROUTINE SAVE 
C 

c *** This routine saves input data on permanent file tape2 
c 

C0MM0N/MCSm/C0RP<78) .ANTM<2.12) «V(25) »NMEAS«KN* 

1 SOROOT (78) f X (25 ) » IU*NTmEAS» IC0UNT» 

2 INAME (<») < ITMEAS»NIMEAS»NTEMP 
COMMON/MISC/IOINf lOOUT, ICARD 
REWIND NTmeaS 

C ***SAVE NOW OF MEAS.« ELEMENTS IN COkR MATRIX. SEX FLAG 
WRITE (NTMEAS) NMEAS.KN. INAME 
C *** SAVE MEAN/STD. vector. 

WRITE (NTMEAS) ( ( ANTM ( 10. JO) . 10=1 .2) . JO= 1 .NmeAS) 

C *** SAVE CORRELATION MATRIX 

WRITE (NTMEAS) (CORR ( 10) . 10=1 .KN) 

REWIND NTMEAS 
WRITE (6.100) 
return 



C 

100 



C 

c 

c 

c 

c 

c 

c 

c 

c 



format (" MEASUREMENT DATA SAVED ON FILE UM£aS OR TAPE2") 
END 

SUBROUTINE SIMSUB (MTYPE) 

*** FOR A GIVEN SOUAR ROOT MATRIX R OF ORDER M IN UPPER 
*** TRIANGULAR FORM OF A CORRELATION OR A COVARIANCE 

*** Matrix of m columns, this routine computes a vector 

*** V OF PSEUOO SCORES OF LENGTH M SUCH THAT IT WILL 
•** BE A OUASI-RANOOM SAMPLE FROM A NORMAL POPULATION 
«»« WITH A MEAN VECTOR OF ZEROS AnD ThE SAME CORRELATION 
•** OR covariance matrix. I.E. from an N(O.R) MAIRIX. 



COMMON/MCSM/CORR (7ft) .ANTM(P. 1?) .V(?5) .NMFAS.KN. 
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1 S«ROOT(70) ,X(2S) tlU.NTMEAS.lCOUNT. 

2 INAME (A) ,ITMEaS.NlMEAS.NTEMP 
COMMON/M I SC/ lOIN » lOOUT, ICA«0 
OIMENSION Y(A5) 

IJ=0 

00 100 I=1«NMEAS 
V(I)=0.0 

CALL RAN02 (RNORM.MTYPEtOUANtlJ) 

100 Y(I)=PNORM 

00 110 I=1*NMEAS 
00 no J=1.I 

K=NOX( J,I) 

no v(I)=v{i)»y(j)*soroot(k) 

RETURN 
END 

SUBROUTINE SORMX (NERROR) 

FOR A GIVEN SYMMETRIC MATRIX OF ORDER N STORED IN ITS 

*** UPPER triangular eohm» columnwise in the vector w or 
*** length N*(N*1)/2. this routine transforms it into its 

"RQUARE ROOT" IN THE SAME FORM IN THE '/ECTOR S0R00T(78) 

COMMON/mCSM/CORR (7fl) .ANTM(2«12) »V(25 ) »Nm£aS.KN» 

1 SQROOT(78) .X (25) . lU .NTME AS » I COUNT . 

2 INamE( 4) ♦iTMEAStNiMFAStNTEMP 
C0MM0N/MISC/I0 IN» ioout, icaro 
DIMENSION ZX(7) 

INTEGER ZX 

C ••*NFRR0P flag is used to SIGNIFY A CONDITION 

C ***WHICH WILL produce a negative SOROOT(N) value 
c ***JUST SEF'OHE implementation OF ThE SQUARE ROOT 
C ***FUNCTION. 

NFRR0R=0 

00 120 I=1,NMEAS 
ZX (4)=I 
N=NOX(I,I) 

C ***TFST FOR negative RADICAL ARGUMENT 
IF (SOROOT (N) .LT.O) GO TO 130 
SO=SORT (SOROOT (N) ) 

DO 100 J=I,NMEAS 
N=NOX(T,J) 

SOROOT <N) =SOROOT (N) /SO 
100 CONTINUE 

C ***NORMAL FXIT FROM SUBROUTINE 
IF (I.FO.NMFAS) NERRORsO 
IF a.FO.NMEAS) RETURN 
I1 = I*1 
ZX(5)=T1 

DO no I2=I1*NmEaS 
00 no J2=I2»NMEAS 
N=NDX { I2*J2) 

L=NDX(I.I2) 

K=N0X(I,J2) 

SOROOT (M) =SOROOT (N) -SOROOT (L) *S0R00T (K) 

ZX(6)=I2 
ZX(7)=J2 
no CONTINUE 
ZX(1)=N 
ZX(2)=L 
ZX(3)=K 

120 CONTINUE 
130 CONTINUE 

c ***apnormal Exit from subroutine 
nfrroh=i 
WRITE (G.180) 

140 rONTTNUF 
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WPITE (/S»190) 

CAUL INPT (4,1, NVALSZ. IANS. XVAL.IEOF) 

GO TO (170, 160.150. 170. 160, 170). 1E0F*1 
150 CONTINUE 

WPITE (6.200) 

WRITE (6»210) (ZX(I) .1=1.7) ,S0R00T(N> 

160 CONTINUE 
RETURN 

170 CC)NT1NUE 

WRITE (6.220) 

GO TO 140 
C 

180 format (•• ANTHROPOMETRIC CORRELATION VALUES ARE INVALID") 

190 format {" — ENTER CARRIAGE RETURN TO CONTINUE PROCESSING"./ 

♦ " — enter ? TO DISPLAY SORMX VARIABLES") 

200 format (" AN INDEFINITE QUANTITY HAS bEEt^l DETECTED BY SORMX"/ 

♦ " variable list within SORMX:"/" N L". 

♦" K I II 12 J2 SOROOT"./.10X,44(1H-) ) 

210 format (8X,4I5.1X.3I5.F9.3,/,10X.44(1H-),//) 

220 format (/" invalid response **♦"/) 

END 

FUNCTION NDXd ,J) 

IF (I.GT.J) go to 100 

NDX=I»(J*J-J)/2 

RETURN 

100 NDX=J»(I*I-I)/2 

RETURN 
END 

FUNCTION UNEORM (mTyPE) 

•••UNIFORM Random number generator on interval (O.d 

COMMON /mCSM/ C0RR(78) ,ANTM(2,12) .V(25) ,NMEAS*KNt 
1S0R00T(7B) ,X(2 S) . IU.NTmEaS. ICOUNT. 

2INAMF(a.) .ITmeaS.NImeaS.NTEMP 
COMMON/MISC/IOIN.IOOUT, ICARD 
c 60 BIT UNIFORM RANDOM NUMBER GENERATOR 

DATA I X/1 677721 3/. C/28 14749767 10655./ 

IF (MTYPE. EO, 32) GO TO 100 

IU=IU*IX 

UNFORM=iU/C 

RFTUPN 

C 32 BIT UNIFORM RANDOM NUMBER GENERATOR 

100 IU=IU^65539 

IF (iu.ge.o) go to no 

IU=1U*2147483647*1 
110 UNFORMsiU 

UNFORM=UNFORM^.4656613E-09 

return 

END 

SUBROUTINE HELP (I) 

RETURN 

END# 
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APPENDIX D 



Program Listing for CAPE Model 
Modified for Use with FORTRAN IV* 



COMMON/MCSM/CORR (78) ,fiNTM(2,12) ,V (25) .NMEfiS.KN, 

1 SQROOT (78) ,X(25) , lU, NTHERS, 1 COUNT, 

2 INfiME (10) ,lTMEflS,NlMEflS,MflXSflM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICflRD 

DflTfi NUCREW/8/ 

DflTfi MflXFUN/4/ 

DflTfi MflX0PT,MflXflNT,MflXflCT/3,7,3/ 

100UT*0 
lSflMP=0 
MMEflS = 12 
100 CONTINUE 
I01N=5 

WRITE (6,290) 

CALL INPT (1, l,NVflLSZ,NMEflS,XVflL, lEOF) 
IEOFPl-lEOF+1 

GOTO (130,100,110,120,280,120) , lEOFPl 
110 CONTINUE 

CALL HELP (1) 

GO TO 100 
120 CONTINUE 

WRITE (6,340) 

GO TO 100 
130 CONTINUE 

IF (NMEflS.LT.l. OR. NMEflS.GT. MMEflS) GO TO 120 
KN*NMEflSM (NMEflS+1) /2 
IINPT^O 
140 CONTINUE 
101N=5 

150 CONTINUE 

WRITE (6,300) 

CALL INPT (l,l,NVflLSZ,IOPT,XVflL,IEOF) 
IE0FP1=IE0F+1 

GO TO (180,140,160,170,280,170) ,IE0FP1 
160 CONTINUE 

CALL HELP (2) 

GO TO 150 
170 CONTINUE 

WRITE (6,340) 

GO TO 150 
180 CONTINUE 

IF (lOPT.LT.l.OR.IOPT.GT.MflXflNT)GO TO 170 
GO TO (190,200,220,240,260,280,225) ,10PT 
190 CONTINUE 

CALL flCCPT (I SAMP) 

IINPT=1 
GO TO 140 
200 CONTINUE 

IF(IINPT.EQ.l) GO TO 210 
WRITE (6,350) 

GO-TO 140 
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210 CONTINUE 
CALL EDM 
GO TO 140 
220 CONTINUE 

IF (ISAMP.EQ. 1) GO TO 230 
222 WRITE (6,320) 

GO TO 140 
225 CONTINUE 

IF (ISflMP.NE. 1) GO TO 222 
CALL PUNRT 
GO TO 140 
230 CONTINUE 
CALL SAVE 
GO TO 140 
240 CONTINUE 

IF (IINPT.EQ. 1) GO TO 250 
WRITE (6,360) 

GO TO 140 
250 CONTINUE 
CALL PRINT 
GO TO 140 
260 CONTINUE 

IF (ISAMP.NE.l) GO TO 270 
CALL OPSAMP (ISAMP) 

GO TO 140 
270 CONTINUE 

WRITE (6,330) 

GO TO 140 
280 CONTINUE 

WRITE (6,310) 

STOP 

290 FORMAT (' ENTER NUMBER OF MEASUREMENTS (1-12)2— *) 

300 FORMAT (’ OPT I ONS ( 1 - I NPUT ; 2-EOI T ; 3-SAVE ; 4-PRI NT ; 5-GENERATE ; • 

1 , *6-END;7-PUNCH) •) 

310 FORMAT (5X, 'END OF OPERATOR SAMPLE') 

320 FORMAT (IX, '«k«mCSM140 UNABLE TO SAVE MEASUREMENT OATA - FATAL 
1 ERRORS') 

330 FORMAT (lOX, '«x«MCSM165 UNABLE TO GENERATE SAMPLE - FATAL ERRORS') 
340 FORMAT (/, 'x«xMCSM024 INVAL 10 RESPONSE ') 

350 FORMAT ( 1 OX, ' «««MCSM1 30 ANTHROPOMETRIC DATA MUST BE INPUT BEFORE', 
1/,13X,'THE EDIT FUNCTION.') 

360 FORMAT (1 OX, '«x«MCSMl50 ANTHROPOMETRIC OATA MUST BE INPUT BEFORE', 
1/,13X,’THE PRINT FUNCTION.') 

END 
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SUBROUTINE INPT (I TYPE » I VflLSZ . NVflLSZ, I VflL . XVflL, lEOF) 

COHMON/HISC/IOIN, lOOUT, ICARD 

DIMENSION ICHAR (80) , IVAL (IVALSZ) ,XVAL (IVALSZ) 



DIMENSION NVAL (20) 

LOGICAL M 1 NCHAR(80) 

EQUIVALENCE (NVAL (1) » NCHAR (1)1 
C 
C 

DATR IBLK, IQUES, IDOL, IQUO, IPOST, IPD, ICOM, IMIN, IPLUS/IH 

1 1H$, IHMHj , IH. , IH. , \ 

DATR mBLK/4H / 

ISTOP^O 

IF (IVALSZ. GE.O) GO TO 100 
ISTDP»1 

IVAISZ=IABS (IVALSZ) 

100 CONTINUE 
10UOTE*0 
NVALSZ*0 
IEOF-0 

IF (lOOUT.EQ.n GO TO 120 
DO no I»l,80 
ICHRR (I) =IBLK 
no CONTINUE 

READ (lOIN, 330, END=ni) ICHAR, NCHAR 
DO 109 J»l,80 

IF (ICHAR (J) ,EQ. IBLK) GO TO 109 
GO rO 130 
109 CONFINUE 
GO rO 113 
in CONFINUE 

REHINO lOlN 
113 lEOf = 1 
RETURN 

120 CONFINUE 
I1«ICAR0 
GO rO 140 
130 CONFINUE 
1 1 =1 

140 CONTINUE 
I2=B0 

150 CONTINUE 

C HRITE (6,501) lOEC, ISTOP, IVALSZ, NVALSZ, 11,12, lEOF 
IF (ISTOP.EQ. 1. AND. IVALSZ. EQ. NVALSZ) RETURN 
IF (I1.LE.I2) GO TO 160 
IF (NVALSZ. LE. IVALSZ) RETURN 
IE0f*5 
RETURN 

160 CONTINUE 

DO 170 1=11,12 
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IF nCHflR(I) .EQ. IBLK ) GO TO 170 
1F1RST=I 
GO TO 180 
170 CONTINUE 

IF (NVflLSZ.EQ.O) lEOF-1 
IF (NVfiLSZ.GT. IVflLSZ) 1E0F=5 
RETURN 

180 CONTINUE 

IF (ICHfiR(IFIRST) .EQ.IQUES) IE0F=2 
IF UCHflR(IFIRST) .EQ.IDOL) IE0F=4 
IF (lEOF.NE.O) RETURN 
IF (ITYPE.EQ.4) GO TO 300 
IF (ITYPE.NE.3) GO TO 190 

IF (ICHflR(IFIRST) .NE.IQUO.RND.ICHflR(IFIRST) .NE. IPOST) GO TO 190 
IFIRST=IFIRST+1 
IQU0TE=1 
190 CONTINUE 

DO 220 I*IFIRST,I2 
IF (IQUOTE.EQ.O) GO TO 200 

IF (ICHflR(I) .EQ.IQUO.OR.ICHfiR(I) .EQ.IPOSTIGO TO 210 
GO TO 220 
200 CONTINUE 

IF (ICHflR(I) .NE.IBLK.flND.ICHfiR(I) .NE.ICOM) GO TO 220 
210 CONTINUE 
lLflST-I-1 
GO TO 230 
220 CONTINUE 
lLflST=I2 
230 CONTINUE 

NVflLSZ=NVflLSZ+l 

Il»ILflST^2 

IF (NVfiLSZ.GT. IVRLSZ) GO TO 150 
IF (ITYPE .EQ.3) GO TO 310 
ISIGN-1 

IF (ICHfiRdFlRSTJ .NE.lMlN.fiND.ICHRR(IFlRST) .NE.IPLUS) GO TO 240 
IF (ICHfiRUFlRST) .EQ.IMIN) IS1GN=-1 
IFIRST=IF1RST+1 
240 CONTINUE 
1NUM»0 
10EC=0 
10ECF=0 

00 260 I=1F1RST, ILfiST 
IF (ICHfiR(l) .NE.IPD) GO TO 250 
I0ECF«IDECF*^1 
IOEC=I 
GO TO 260 
250 CONTINUE 

CfiLL ICONV (ICHRR m . IfiDD, lEOF) 

IF (IE0F.EQ.3) RETURN 
lNUM=INUMHlO-^lfiDO 
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501 FORMAT {/ • TEST \ lOIlO) 

IF (lOEC.NE.O) I0EC=IDEC+1 
260 CONTINUE 

IF (IOECF.lt. 2) GO TO 280 
270 CONTINUE 

C HRITE (6,501) lOEC, lOECF 
IE0F*3 
RETURN 

280 CONTINUE 

INUM«INUM«ISIGN 
XNUM=FLOAT (INUM) 

CC WRITE (6,501) lOEC, IDECF, INUM, lAOO, ISIGN, ITTPE 

IF ( (ITTPE.EQ. I.RNO. lOEC.LE. 1) .OR. (ITTPE.EQ.2) ) GO TO 290 

1E0F«3 

RETURN 

290 CONTINUE 

IF (lOEC.NE.O) XNUM=XNUM/10.0mk (lOEC-n 
IF (ITTPE.EQ. 1) IVAL (NVALSZ) =IFIX (XNUM) 

IF (ITTPE.EQ. 2) XVRL (NVALSZ) =XNUM 
ICRR0=ILRST-^1 

C WRITE(6,50n I0EC,I0ECF,INUM. IRDO, I S 1 GN, IT YPE, NVALSZ, 

C 1 IVAL (NVALSZ) , ICRRO 
GO TO 150 
300 CONTINUE 

ILAST=IVALSZh4 

909 FORMAT! 2X, 80A1) 

910 FORMAT (2X, 20A4) 

310 CONTINUE 

DO 320 J»1,IVALSZ 
IVAL (J) -I4BLK 
320 CONTINUE 

L=ILAST-IFIRST+1 
NVALSZ* (L-1) /lO+l 
IF (L.GT.4 hIVRLSZ) go TO 270 
C ENCOOE (L, 340, IVAL (1) ) (ICHAR (I) , IFIRST, ILAST) 

J - 0 

904 FORMAT (• ENCOOE \ 101 10) 

00 311 I » IFIRST, ILAST . 

J - J ♦ 1 

NCHAR(J) = NCHAR(I) 

311 CONTINUE 

DO 313 I * 1,IVALSZ 
IVAL (I) *NVAL(I) 

313 CONTINUE 

1 CARO* II 
RETURN 

330 FORMAT (80A1, Tl, 80A1) 

340 FORMAT (80A1) 

END 
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SUBROUTINE ICONV t JCHflR, I M EOF) 

DIMENSION JNUM(IO) 

OflTfl JNUM/lHO, IHl , 1H2, 1H3» IHM, 1H5* 1H6, 1H7* 1H8, 1H9/ 
11 FORMAT (//• RT. ICONV \ A2,2I10 /) 

DO 100 J»l. 10 

IF (JCHAR.NE. JNUM (J)l GO TO 100 

I1*J-1 

RETURN 

100 CONTINUE 
lEOF - 3 

WRITE (6,11) JCHAR, II, lEOF 

RETURN 

END 
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SUBROUTINE flCCPT (ISflMP) 

COMMON/MCSM/CORR (78) ,flNTM(2,l2) , V (25) .NMEflS.KN, 

1 SQROOT (78) ,X (25) , IU,NTMEfiS, I COUNT. 

2 I NOME (10) . ITMEflS.NiMEfiS.MflXSflM.NTEMP 
COMMON/MISC/IOIN. lOOUT. ICflRO 
DIMENSION BUFF (3) .TEMP (12) 

ORTfl BUFF/3 hO./ 

ISflMP=0 

IC0UNT*0 

INL1NE=0 

10IN*5 

100 CONTINUE 

WRITE (6,610) 

CALL INPT (1. l.NOINPT. IMOOE.XINPT, lEOF) 
IE0FP1=IE0F+1 

GO TO (130. 1 10. 120. 1 10.490. 1 10) . lEOFPl 
no WRITE (6.650) 

GO TO 100 
120 CALL HELP (4) 

GO TO 100 

130 IF (IMOOE.EQ. 1) GO TO 200 
IF (IM00E.NE.2) GO TO 110 
140 CONTINUE 

WRITE (6.620) 

CALL INPT (1, l.NOINPT. IFFORM.XINPT, lEGF) 
IE0FP1=IE0F+1 

GOTO (170,150,160,150.490,150). lEOFPl 
150 WRITE (6,650) 

GO TO 140 
160 CALL HELP (5) 

GO TO 140 

170 IF (IFFORM.EQ. 1) GO TO 190 
IF (IFF0RM.NE.2) GO TO 150 
REWIND NTMEAS 

READ (NTMEAS. END=180) NM^ AS, KN , I NAME 

READ (NTMEAS) ( (ANTM (10. JO) , 10=1 , 2) , J0=1 , NMEAS) 

READ (NTMEAS) (CORR ( 1 0) , 1 0= 1 , KN) 

REWIND NTMEAS 
WRITE (6.690) 

ISAMP=1 

RETURN 

180 CONTINUE 

WRITE (6.670) 

GO TO 490 
190 I0IN=1 
200 CONTINUE 

CALL NAMEFL (lERROR) 

ICOUNT*ICOUNT+IERROR 
C WRITE (6,901) IMOOE. lERROR, ICOUNT 
901 FORMAT!/ ’ ACCPT \ lOIlO) 
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INLINE*INLINE^1 
IF (IMODE.EQ.l) WRITE (6,500) 

IF (IMOOE.NE.n WRITE (6,510) INflME 

MflXlNP*2 

C WRITE(6.90!) IMODE, lERROR, ICOUNT, INLINE. MflXiNP, lEOF 

DO 320 IROW=l.NMEflS 
210 CONTINUE 

IF (IMOOE.EQ.l) WRITE (6,530) IROW 

220 NflMT=0 
230 CONTINUE 

lTOTlN=MflXlNP-NflMT 

IC0L=NRMT+1 

CALL INPT (2. IT0TIN,N0INPT, IVflL, flNTM (I COL, IROW) . I EOF) 

IF (IE0F.EQ.5.0R. lEOF.EQ.O) NflMT=NflMT>NOlNPT 
C WRITE(6,90n IMODE, ITOTIN, MflXINP, NRMT, lEOF, ICOL, ICOL 

IF (IMOOE.EQ.l) GO TO 240 

INLINE=INLINE*^1 
IF (NfiMT.EQ. MflXiNP) GO TO 290 
IF (lEOF.EQ.l) GO TO 480 
WRITE (6,540) IROW 
IF (IE0F.EQ.5) GO TO 270 
WRITE (6,700) INLINE 
IC0UNT»IC0UNT>1 
GO TO 320 
240 CONTINUE 

IE0FP1 = IE0F-^1 

GOTO (280,210,250,260,490,270), lEOFPl 
250 CONTINUE 

CALL HELP (6) 

GO TO 210 
260 CONTINUE 

WRITE (6,650) 

GO TO 210 
270 CONTINUE 

I0LI5T*N0INPT-IT0TIN 
WRITE (6,600) lOLIST 
IF (IMOOE.NE.n IC 0 UNT=IC 0 UNT -^1 
GO TO 290 
280 CONTINUE 

IF (NflMT.GE.2) GO TO 290 
WRITE (6,550) 

GO TO 230 
290 CONTINUE 

IF (NflMT.GT.2) NflMT=2 
IRERRsO 

DO 300 ICOL = UNflMT 

CALL REASON (1 , ICOL. IROW, ANTM (1 , IRON) , ANTM (2, IROW) , IFATAL) 
IRERR=»IRERR+IFATAL 
300 CONTINUE 

IF (IRERR.EQ.O) GO TO 320 
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IF (IMODE.NE. 1) GO TO 310 
NRITE (6,630) IROW 
GO TO 220 

310 ICOUNT = ICOUNT-^IRERR 
320 CONTINUE 

IF (iMODE.EQ.n WRITE (6,560) 

IF (IMOOE.NE.I) WRITE (6,520) 

DO 460 IROW=l,NMEflS 
330 CONTINUE 
NflMT=0 

MflXlNP=NMEfiS+l-IROW 
340 continue 

IF (IMOOE.EQ. l.fiND. IROW.EQ.NMEflS) WRITE (6.580) IROW. MAX I NP. NMEflS 
IF (IMOOE.EQ. 1. AND. IROW. NE.NMEfiS) WRITE (6.570) IROW, MAXINP. IROW, 
INMEAS 

350 CONTINUE 

IT0TIN=MAXINP-NAMT 
ISTART = IROW-»NAMT 

CALL INPT (2,IT0TIN,N0INPT. IVAL,TEMP (ISTART) .lEOF) 

IF (lEOF.EQ.O.OR. IE0F.EQ.5) NAMT = NAMT-^NOINPT 

IF (IMODE.EQ.n GO TO 360 

INLINE«INLINE^1 

IF (NAMT.EQ. MAXINP) GO TO 410 

IF (lEOF.EQ.O) GO TO 350 

IF (lEOF.EQ.l) WRITE (6,710) 

IF (lEOF.EQ. 1) GO TO 470 

IF (IE0F.EQ.5) GO TO 390 

IC0UNT=IC0UNT+1 
WRITE (6,700) INLINE 
GO TO 460 
360 CONTINUE 

lEOFPmlEOF^l 

GO TO (400,380,370,380,490.390) . lEOFPl 
370 CONTINUE 

CALL HELP (7) 

GO TO 340 
380 continue 

WRITE (6,650) 

GO TO 340 
390 CONTINUE 

I0LIST=N0INPT-IT0TIN 
WRITE (6,600) lOLIST 
IF (IMOOE.NE.I) IC0UNT=IC0UNT+1 
NAMT=MAXINP 
GO TO 410 
400 CONTINUE 

IF (NOINPT.EQ. ITOTIN) GO TO 410 
IT0T1N*MAXINP-NAMT 
WRITE (6,590) ITOTIN, IROW 
GO TO 350 
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CONTINUE 

lLfiST=IROH'^NfiMT-l 

1RERR=0 

00 420 ICOL=IROW, ILflST 

WRITE (6,902) IROW, ICOL, ILPST, IMODE, TEMP(ICOL) 

FORMRTI’ fiCC’, 4110, 4F10.3) 

CRLL RERSON (2, 1 COL , 1 ROW, TEMP (I COL) , 0, I FRTRL ) 

IRERR=1RERR+IFRTRL 

WRITE (6,901) IROW, ICOL, IRERR, I FRTRL 
CONTINUE 

IF (IRERR. EQ.O) GO TO 440 
IF (IMODE. EQ. 2) GO TO 430 
WRITE (6,640) IROW 
GO TO 330 

ICOUNT=ICOUNT+IRERR 
GO TO 460 
CONTINUE 

DO 450 1C0L=IR0W,NMERS 
N-NDX (IROW, ICOL) 

CORR (N) -TEMP (ICOL) 

CONTINUE 

CONTINUE 

CONTINUE 

IF (ICOUNT.EQ.O) ISflMP=I 
IF (IMOOE.NE.I) WRITE (6,680) ICOUNT 
GO TO 490 
WRITE (6,660) 

RETURN 

CONTINUE 

RETURN 

FORMRT (5X. 'ENTER PRIRS OF MERNS RNO STRNORRD OEV I RT I ONS-- ') 

FORMRT (//, ' DIAGNOSTICS OF FORMATTED RERO OPER OESC: ',10R4 

1 ./) 

FORMRT (//, IX, 'CORRELRTION/RNTHROPOMETRIC MEfiSUREMENTS: ') 

FORMRT (' MEASUREMENT - ',12,' ?— -') 

FORMRT (' MEASUREMENT -',12) 

FORMRT (5X, 'ENTER STANDARD DEVIATION --') 

FORMRT (5X, 'ENTER ANTHROPOMETRIC CORRELATION MATRIX ORTA IN R ROWW 
USE UPPER ',/,5X, 'TRIANGULAR FORM. ',/, 5X, 'NOTE : ALL DIAGONAL ELEM 
2ENTS MUST EQUAL 1.000') 

FORMAT (5X,4HR0W ,12, '(',12,' ELEMENTS, COLUMNS ' , 12, , 12, ') ’) 
FORMAT (5X,4HR0W ,12, '(',12,' ELEMENT, COLUMN ',12,')') 

FORMAT (10X,23H«hhENTER THE REMAINING ,I2,18H ELEMENTS FOR ROW ,12 
1 ) 

FORMAT (lOX, ' h««ACCPT051 THE LAST ',12,' INPUTS HAVE BEEN IGNORED' 
1 ) 

FORMAT (5X, 'INPUT MODE (1* I NTERACT I VE ; 2=F I LE) — ') 

FORMAT (5X,'REA0 MODE (1=FORMATTEO;2=UNFORMATTED) -- ') \ 

FORMAT (lOX, 'H«HfiCCPT060 RE-ENTER ORTA FOR ',12) 

FORMAT (lOX, 'HHMRCCPT1300 RE-ENTER DATA FOR ROW ',13) 
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650 FORMAT (1 OX. *HMNflCCPT 155 INVALID RESPONSE') 

660 FORMAT (1 OX. 'hhhACCPT UNEXPECTED EOF ON NIMEAS FOR MEAS '.12) 

670 FORMAT (1 OX. 'hh^READFIMOSO FILE ANTSTA IS NOT ATTACHED TO THE PROG 
IRAM'./. lOX. ' hhhUNFORMATTED READ ON THIS FILE IS ABORTED --'//) 

680 FORMAT dOX. ' hhhACCPT3000 TOTAL ERRORS ON INPUT DATA ’.15) 

690 FORMAT (' MEAN. STD DEV 4;C0RRELAT I ON DATA READ') 

700 FORMAT (1 OX. 'mhmACCPT 150 INVALID DATA AT LINE '.15) 

710 FORMAT (1 OX. 'hhxACCPT1200 UNEXPECTED END OF FILE') 

END 
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SUBROUTINE EDIT 

COMMON/MCSM/CORR (78) ,RNTM (2, 12) ,V (25) .NMEfIS, KN, 

1 SQROOT (78) ,X (25) , IU,NT)iEfiS, ICOUNT, 

2 INRME (10) • ITMEflS,NlHEflS,MfiXSflM,NTE)iP 
C0MM0N/MI3C/I0IN, lOOUT, ICflRD 
DIMENSION BUFF (3) 

DRTfl BUFF/3 hO./ 

I0IN«5 

100 CONTINUE 

WRITE (6,350) 

CALL INPT (l,l,NVALSZ,IOPT,XVflL,IEOF) 
1E0FP1=IE0F>1 

GOTO (130,100,110,120,340,120), lEOFPl 
no CONTINUE 

CALL HELP (8) 

GO TO 100 
120 CONTINUE 

WRITE (6,360) 

GO TO 100 
130 CONTINUE 

IF (lOPT.LT. l.OR. I0PT.GT.5) GO TO 120 
GOTO (140, 140,240,330,340) . lOPT 
140 CONTINUE 
ICOL-IOPT 
150 CONTINUE 
MAXINP*2 
NAMT*0 

160 CONTINUE 

IF (ICOL.EQ.l) WRITE (6,380) 

IF (IC0L.EQ.2) WRITE (6,390) 

170 CONTINUE 

IVAL3Z==MAXINP-NAMT 

JsNAMT*^! 

CALL INPT (2, IVALSZ,NVALSZ, IVAL,BUFF (J) , lEOF) 
IE0FP1»IE0FP1 

GOTO (220,180,200,190,100,210), lEOFPl 
180 CONTINUE 

WRITE (6,370) 

GO TO 160 
190 CONTINUE 

WRITE (6,360) 

GO TO 160 
200 CONTINUE 

IF (ICOL.EQ.l) CALL HELP (11) 

IF (IC0L.E0.2) CALL HELP (12) 

GO TO 160 
210 CONTINUE 

lOLIST^NVALSZ-IVALSZ 
WRITE (6,400) lOLIST 
220 CONTINUE 
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IROW=BUFF m 
NftMT=NVRLSZ+NRMT 
IF (NRMT.GE.MRXINP) GO TO 230 
IF nCOL.EQ.n WRITE (6,420) IROW 
IF (IC0L.EQ.2) WRITE (6,430) IROW 
GO TO 170 
230 CONTINUE 

CRLL PERSON (1,IC0L,IR0W,BUFF (2) ,0,IFRTRL) 

IF (IFRTRL.EQ.l) GO TO 150 
RNTM (ICOL, IROW) =BUFF (2) 

GO TO 100 
240 CONTINUE 
MRXINP-3 
NRMT*0 

250 CONTINUE 

WRITE (6,410) 

260 CONTINUE 

IVRLSZ*MRXINP+NRMT 

J»NRMT+1 

CRLL INPT (2, IVRLSZ,NVRLSZ. IVRL,BUFF (J) , lEOF) 
IE0FP1*IE0F+1 

GOTO (310,270,290,280, 100,300) , lEOFPl 
270 CONTINUE 

WRITE (6,370) 

GO TO 250 
280 CONTINUE 

WRITE (6,360) 

GO TO 250 
290 CONTINUE 

CALL HELP (13) 

GO TO 250 
300 CONTINUE 

IOLIST*NVRLSZ-IVRLSZ 
WRITE (6,400) lOLIST 
GO TO 320 
310 CONTINUE 

ICOL*BUFF (1) 

IROW=BUFF (2) 

NRMT^NRMT-^NVRLSZ 
IF (NRMT.GE.MRXINP) GO TO 320 
IF (NRMT.EQ.l) WRITE (6,440) ICOL 
IF (NRMT.EQ.2) WRITE (6,450) ICOL, IROW 
GO TO 260 
320 CONTINUE 

CRLL RER50N (3, 1 COL, I ROW. BUFF (3) , 0. IFRTRL) 

IF (IFRTRL.EQ.l) GO TO 240 

CRLL PERSON (2, I COL , I ROW, BUFF (3) , 0, IFRTRL) 

IF (IFRTRL.EQ.l) GO TO 240 
N=NOX (IROW, ICOL) 

CORR (N) =BUFF (3) 
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GO TO 100 
330 CONTINUE 
COLL PRINT 
GO TO 100 
3M0 CONTINUE 
RETURN 

350 FORMAT (SX.'EDIT MODE (1 -MEANS; 2-STO ; 3-flNTH. CORR; M-PR I NT ; 5-END) — M 
360 FORMAT ( 1 OX , ’mmmEO I T04M INVALID RESPONSE') 

370 FORMAT (1 OX, 'mxmEO I T095 NO INFORMATION ENCOUNTERED’) 

380 FORMAT (5X, ’INPUT INDEX AND MEAN VALUE’) 

390 FORMAT (5X, ’INPUT INDEX AND STANDARD DEVIATION’) 

^00 FORMAT (lOX, ’ mn«EOIT 0M2 THE LAST’,12,’ INPUTS HAVE BEEN IGNORED’) 
mo FORMAT (5X, ’INPUT COLUMN , ROW AND CORRELATION’) 

420 FORMAT (5X, ’INPUT MEAN VALUE FOR INDEX ’,13) 

430 FORMAT (5X, ’INPUT STANDARD DEVIATION FOR INDEX ’,13) 

440 FORMAT (5X, ’INPUT ROW AND CORRELATION FOR COLUMN M2) 

450 FORMAT (5X, ’INPUT CORRELATION FOR ROW ’,12,’ COLUMN ’,12) 

END 
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SUBROUTINE NftMEFL (lERRORI 

COMMON/MCSM/CORR (78) .RNTM (2, 12) , V (25) .NMEPS.KN, 

1 SQROOT (78) ,X (25) , lU.NTMEflS, ICOUNT, 

2 INPME (10) , ITMEfiS.NIMEPS.MflXSflM.NTEMP 
COMMON/MISC/IOIN, lOOUT, ICflRD 
DIMENSION IBUFF (10) ,XVflL (10) 

DflTfi IBUFF/IOmMH / 

DPTfl lyBLK/MH / 
lERRORsO 
DO 100 j*i,y 
IBUFF U) =IMBLK 
100 CONTINUE 
110 CONTINUE 

IF (I0IN.E0.5) WRITE (6,180) 

CRLL INPT (M, 10, NVRLSZ, IBUFF, XVfiL. lEOF) 

IF (IE0F.LT.2.0R. IE0F.EQ.5) GO TO 150 
IF (I0IN.NE.5) GO TO IMO 
IEOFMUIEOF-1 

GO TO (120,130, 170) , lEOFMl 
120 CONTINUE 

COLL HELP (ly) 

GO TO no 
130 CONTINUE 

WRITE (6,190) 

GO TO no 
lyo CONTINUE 

WRITE (6,190) 

IERROR=I 
150 CONTINUE 

DO 160 J>1, 10 
INRME (J) =IBUFF (J) 

160 CONTINUE 
170 RETURN 

180 FORMRT (5X, 'OPERATOR SAMPLE DESCRI PTI ON (MO CHAR MAX) ’) 

190 FORMAT (1 OX, *hhkNAMEFL027 I NVAL I D RESPONSE ' ) 

END 
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SUBROUTINE OPSflMP (ISflMP) 

COMMON/MCSM/CORR (78) ,RNTM (2, 12) .V (25) .NMERS.KN, 

1 SQROOT (78) ,X (25) • lU.NTMERS. I COUNT, 

2 INRME (10) , ITMERS.NIMERS.MRXSRM.NTEMP 
COMMON/MISC/IOIN, lOOUT, ICRRD 
DIMENSION IRRN(9) 

DRTR IRRN (1) , IRRN (2) , IRRN (3) , IRRN (4) , IRRN (5) , IRRN(6) , IRRN (7) , 
lIRRN (8) , IRRN (9) Z33333333, 55555555, 77777777, 222222221 , 

1 444444443,666666669.888888889,999999999. 123456789/ 

I0IN«5 

DO 100 K==1,KN 
SQROOT (K) =C(3RR (K) 

100 CONTINUE 

CRLL SQRMX (NERROR) 

IF (NERROR. EQ. 1) RETURN 
no CONTINUE 

WRITE (6,320) 

CRLL INPT (1, 1,NVRLSZ.NSRMP,XVRL, lEOF) 
lEOFPl-IEOF^l 

GO TO (150, 120, 130, 140,290. 140) , lEOFPl 
120 CONTINUE 

WRITE (6,330) 

GO TO 110 
130 CONTINUE 

CRLL HELP (15) 

GO TO 110 
140 CONTINUE 

WRITE (6,340) 

GO TO no 
150 CONTINUE 

IF (NSRMP.GT.O.RND.NSRMP.LE.MRXSRM) GO TO 160 
WRITE (6,350) MRXSRM 
GO TO no 
160 CONTINUE 

WRITE (6,300) 

CRLL INPT (1, 1,NVRLSZ.MTTPE,XVRL, I EOF) 

1E0FP1»IE0F'^1 

GO TO (190,160,170,180.290,180) ,IE0FP1 
170 CONTINUE 

CRLL HELP (20) 

GO TO 160 
180 CONTINUE 

WRITE (6,340) 

GO TO 160 
190 CONTINUE 

IF (MTTPE.NE.32.RND.MTTPE.NE.60) GO TO 180 
200 CONTINUE 

WRITE (6,310) 

CRLL INPT (1, 1,NVRLSZ, IS,XVRL, lEOF) 

IE0FP1=IE0F+1 
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GOTO (240,210.220,230,290,230) , lEOFPl 
210 CONTINUE 

WRITE (6,330) 

GO TO 200 
220 CONTINUE 

CALL HELP (21) 

GO TO 200 
230 CONTINUE 

WRITE (6,340) 

GO TO 200 
240 CONTINUE 

IF (IS.EQ.O) GO TO 250 
IF (IS.LT. l.OR. IS.GT.9) GO TO 230 
IU«IRfiN (IS) 

GO TO 260 
250 CONTINUE 

IU=67111133 
260 CONTINUE 

REWIND NTEMP 

WRITE (NTEMP) INflME.NMEfiS 
DO 280 I=l,NSnMP 
CALL SIMSUB (MTYPE) 

00 270 K=1,NMEAS 
V (K) =V (K) mANTM (2,K) +ANTM (1 ,K) 

CONTINUE 



270 

280 



WRITE (NTEMP) 
CONTINUE 
ISAMP-1 

ENO FILE NTEMP 
REWIND NTEMP 
REWIND NTMEAS 
WRITE (6,370) 
370 FORMAT ( lOX, 



(V (10) , I0=1.NMEAS) 



MH SAMPLE DATA SAVED ON FILE 9M 





WRITE 


(6, 360) 


290 


CONTINUE 

RETURN 


300 


FORMAT 


(5X. ’E 


310 


FORMAT 


(5X, ’E 


320 


FORMAT 


(5X, M 


330 


FORMAT 


(lOX, ’ 


340 


FORMAT 


(lOX, ’ 


350 


FORMAT 


(lOX, ’ 


360 


FORMAT 

ENO 


(lOX, ’ 



(32 OR 60 BIT) — 
IF NONE OR 1-9) — 



MACHINE WORD SIZE 
RANDOM NO SEED (0 
SAMPLE SIZE — M 

^SAMP215 NO INFORMATION ENCOUNTERED*) 
"SAMP411 INVALID RESPONSE’) 

"SAMP009 SAMPLE SIZE MUST BE BETWEEN 1 AND 
^MPLE DATA SAVED ON FILE SAMPLE OR TAPE 3.’) 



. 13 ) 
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SUBROUTINE PRINT 

COMM0N/MCSM/CORR (78) ,PNTM (2, 12) ,V (25) ,NMEftS,KN, 

1 SQROOT (78) ,X (25) , lU.NTMEflS, ICOUNT, 

2 iNflME (10) , lTMEftS,NlMEflS,MftXSPM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICflRD 

INTEGER FMT.FMTl 

DIMENSION FMT (6) ,NUM (12) .TEMP (12) , INDEX (2) ,FMTl (5) 

DPTR FMTl /’(IX, \ •(2X.’, ’R2.2’, ’X))V 

DflTR FMT /MIX, ’ (6X) ’, ’ ’, ’F6.3’,’)’/ 

ORTR NUM/lHl, 1H2, IH3, 1H4, 1H5, 1H6, 1H7, 1H8, 1 H9, 2H10, 2H1 1 , 2H1 2/ 
FMTl (2) *NUM(NMERS) 

IST0RE=0 
I0FLRG=0 
100 CONTINUE 

WRITE (6,630) 

CRLL INPT (1,1,NVRLSZ,I0PT,XVRL,IE0F) 

IE0FP1=IE0F+1 

GO TO (120, 1 10, 1 10, 1 10,390, 1 10) , lEOFPl 
no CONTINUE 

WRITE (6,650) 

GO TO 100 
120 CONTINUE 

IF (10PT.LT.1.0R.I0PT.GT.4)G0 TO 110 

IF (lOPT.EQ.l) GO TO 130 

IF (I0PT.EQ.2) GO TO 150 

IF (10PT,EQ.4) GO TO 390 

GO TO 190 
130 CONTINUE 

WRITE (6,580) 

DO 140 1R0W=1,NMERS 

WRITE (6,640) IROW. (RNTM (JO, IROW) , J0=1 .2) 

140 CONTINUE 

WRITE (6,620) 

GO TO 100 
150 CONTINUE 

WRITE (6,590) 

WRITE (6, FMTl) (NUM (10) , 10=1 ,NMERS) 

WRITE (6,600) 

DO 180 IR0W=1,NMERS 
DO 160 1C0L=1,NMERS 
N=NDX (IROW, ICOL) 

TEMP (ICOL) =CORR (N) 

160 CONTINUE 

IF (IROW.GT.l) GO TO 170 

WRITE (6,610) (TEMP (ICOL) , I COL = I ROW, NMERS) 

GO TO 180 

170 FMT (2) «NUM (lROW-1) 

FMT (4) =NUM (NMERS+l-IROW) 

WRITE (6, FMT) (TEMP (ICOL) , ICOL = IROW, NMERS) 

180 CONTINUE 
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WRITE (6,620) 

GO TO 100 
190 CONTINUE 

WRITE (6,430) 

IF (IOFLRG.lt. 1) WRITE (6.440) 
lOFLRG*! 

CRLL INPT (1,2.NVRLSZ. INDEX.XVRL, lEOF) 
IE0FP1*IE0F-^1 

GOTO (230.270,200.210.100.220). lEOFPl 
200 CONTINUE 

CALL HELP (19) 

10FLRG=0 
GO TO 190 
210 CONTINUE 

WRITE (6,420) 

GO TO 190 
220 CONTINUE 

I0LlST=NVRLSZ-2 
WRITE (6,450) lOLIST 
230 CONTINUE 

IF (INOEX (1) .GT.O) GO TO 240 
WRITE (6,480) 

GO TO 190 
240 CONTINUE 

IF (ISTORE.GT.O.RNO. INOEX (1) .GT. ISTORE) GO TO 250 
NRSKIP=INOEX (1) -1 
GO TO 260 
250 CONTINUE 

NRSKIP=INDEX (1) -ISTORE-1 
260 CONTINUE 

ISTRRT=INOEX (1) 

IF (NVRLSZ.EO. 1) GO TO 280 
IF (NVRLSZ.GE.2) GO TO 290 
270 CONTINUE 
ISTORE-0 
ISTRRT*! 

IEN0*999 
NRSKIP*0 
GO TO 320 
280 CONTINUE 

IENO=ISTRRT 
GO TO 310 
290 CONTINUE 

IF (INOEX (1) .LE. INDEX (2) ) GO TO 300 
WRITE (6,490) 

GO TO 190 
300 CONTINUE 

lENO^lNDEX (2) 

310 CONTINUE 

IF (INOEX (1) ,GT. ISTORE. RNO. ISTORE. NE.O) GO TO 330 
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320 CONTINUE 

REWIND NTEMP 

RERD (NTEMP, EN0=41D) INRME.NMERS 
330 CONTINUE 

IF (NRSKiP.LT.n GO TO 350 
DO 340 I=1,NRSKIP 

RERD (NTEMP, END=4DD) (V (K) , K= 1 , NMERS) 

340 CONTINUE 
350 CONTINUE 

WRITE (6,530) (I NRME ( 10) , 1 0= 1 , 1 0) 

IF (ISTRRT.EQ. lENO) WRITE (6,540) lENO 

IF (ISTfiRT.NE.IEN0.RN0.IEN0.NE.999) WRITE (6.550) ISTRRT.IENO 
WRITE (6.500) 

WRITE (6,FMT1) (NUM (10) , 10=1 .NMEfiS) 

WRITE (6,600) 

ISKIP=0 

NF0UN0=0 

DO 360 I=ISTfiRT, lEND 
IF (ISKIP.EQ.n WRITE (6.460) 

ISKIP*D 

IF (ISTRRT.NE.IENO.RNO.I.NE.IENO.RNO. (I /5) m5. EQ. I ) ISKIP = 1 
RERO (NTEMP, EN0=37D) (V (K) ,K=1 .NMERS) 

NF0UN0=NF0UN0*^1 

ISTORE*! 

WRITE (6,510) (V (K) ,K=1 .NMERS) 

360 CONTINUE 
GO TO 380 
370 CONTINUE 

IF (IEND.NE.999.fiND. ISTRRT.NE. lEND) WRITE (6,470) lEND 
ISTORE»I 

IF (NFOUND.lt. 1) WRITE (6,560) 

380 CONTINUE 

WRITE (6,520) 

GO TO 190 
390 CONTINUE 

REWIND NTEMP 
RETURN 

400 CONTINUE 

I0LIST=I"2-^IST0RE 
WRITE (6,560) 

GO TO 190 
410 CONTINUE 

WRITE (6,570) 

GO TO 100 

420 FORMRT (1 OX, 'mmhPRINT 1 35 INVRLIO RESPONSE’) 

430 FORMRT (5X, 'OUTPUT MODE; --’ ) 

440 FORMRT (lOX, 'FOR RLL OPERRTOR DRTR (ENTER CRRRIRGE RETURN) ',/, lOX, 
I'FOR INDIVIDURL OPERRTOR DRTfi (ENTER OPERRTOR NUMBER) ',/, 1 OX, 'FOR 
2fl SEQUENCE OF OPERRTORS (ENTER RRNGE OF OPERRTOR NUMBERS) ',/, 1 OX, ' 
3EN0 (ENTER $) ') 
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450 FORMAT ( lOX, INI 1 1 9 THE LAST M4. 'INPUTS HAVE BEEN IGNORED’) 

460 FORMAT (IH ) 

470 FORMAT (1 OX, ’«»<*«PRINT21 1 ERROR IM UPPER INOEX --’.14) 

480 FORMAT (1 OX, ’hmhPR INT 1 20 INDEX NO* LESS THAN OR EQUAL ’,/, ’h><nT 

10 ZERO IS NOT VALID’) 

490 FORMAT ( 1 OX, ’hnhPR INT 1 60 FIRST INOEX MUST BE LESS THAN SECOND INOE 
IX’) 

500 FORMAT (27X, ’OPERATOR SAMPLE ’,/, 71 (1 H-) ) 

510 FORMAT (1 1F6. 2, F5. 2) 

520 FORMAT (71 (1H-) ) 

530 FORMAT (IX, 10A4, lOX, A4) 

540 FORMAT (lOX, ’OPERATOR NO’ , 14) 

550 FORMAT (lOX, ’OPERATOR NOS’, 14,’ --’,14) 

560 FORMAT (lOX, ’n«*<PRINT 1500 INDEX NUMBER EXCEEDS RECORDS IN OPERATOR 
1 FILEhhn’) 

570 FORMAT (1 OX, ’nnhPRINT 1 600 FILE SAMPLE IS NOT ATTACHED TO THE PROGR 
1AM’,/, lOX, ’ hhxOPERATOR OUTPUT REQUEST ABORTED ---’,//) 

580 FORMAT (/, 21X, ’MEANS AND STANDARD DEVI AT IONS’ , /, 1X.70 (1H-) , /, 2X, ’M 
1EASUREMENT’,6X, ’MEANS ’,6X, ’STD DEV’,/, IX, 71 (1H-) ) 

590 FORMAT (/, 7X, ’CORRELAT ION MATRIX FOR ANTHROPOMETRIC MEASUREMENTS’, 
1/, IX, 70 (1H-) ) 

600 FORMAT (1X,70(1H-)) 

610 FORMAT (1X,12F6.3) 

620 FORMAT (1X,70(1H-)) 

630 FORMAT (5X, ’PRINT MODE (l-MEAN/STO; 2-CORR MATRX; ’ . ’3-OPER MEAS;4-EN 
ID) — ’) 

640 FORMAT (7X, 12, 5X, 2 (5X , F7. 2) ) 

650 FORMAT (lOX, ’ mnhPRINT 103 INVALID RESPONSE’) 

END 
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SUBROUTINE RflND2 (RNORM, MTTPE, QUfiN, IJ) 
COMMON/MCSM/CORR (78) ,RNTM (2, 12) ,V (25) ,NMEfiS,KN, 

1 SQROOT (78) ,X (25) , lU^NTMEfiS, ICOUNT^ 

2 iNflME (10) , ITMEflS,NlMEflS,MflXSflM,NTEMP 
DIMENSION QUflN(2) 

COMMON/MISC/IOIN. lOOUT. ICRRD 

1J=IJ^1 

IF (IJ.E0.2) GO TO no 
100 CONTINUE 

X1 = 2.0hUNFORM (MTTPE) -1 .0 
Y1=2.0nUNF0RM (MTYPE) -1 .0 
S=X1nX 1+Y1«Y1 
IF (S.GE.1.0) GO TO 100 
S»SQRT (-2.0KfiL0G(S) /S) 

QURN (1) =X1mS 
QURN (2) =Y1«S 
no CONTINUE 

RNORM=QURN (IJ) 

IF (IJ.E0.2) IJ=0 

RETURN 

END 



104 



0 




SUBROUTINE REASON (ICALL. ICOL. IRON, TSTORK TST0R2, 


1 




COMMON/MCSM/CORR (78) ,RNTM(2, 12) ,V (25) ,NMEfiS,KN, 


2 




1 SQROOT (78) ,X (25) , IU,NTMERS, ICOUNT, 


3 




2 INRME (10) . ITMERS,NIMERS,MRXSRM,NTEMP 


« • 




COMMON/MISC/IOIN, lOOUT, ICRRO 


5 




IFRTRL=0 


6 




1RFLRG=0 


7 


C 


WRITE (6,901) ICRLL, ICOL, IROW, IFRTRL, TSTORl .TST0R2 


8 


901 


FORMAT r REASON \ 4110, 4F10.3) 


9 




IF (ICALt.LT.O) IAFLAG=1 


10 




ICRLL^IRBS (ICRLL) 


11 




IF (ICRLL.GT.I) GO TO 110 


12 




IHRITE»1 


13 




IF (IROW.LT. 1 .OR. IROH.GT.NMEAS) GO TO 120 


14 




IF (IC0L.EQ.2) GO TO 100 


15 




RETURN 


16 


100 


CONTINUE 


17 




IWRlTE-2 


18 




IF (TST0R2.LT.0.) GO TO 120 


19 




lWRlTE-3 


20 




IF (TST0R1.LE.TST0R2) GO TO 120 


21 




RETURN 


22 


1 10 


CONTINUE 


23 




IWRITE-1 


24 




IF (IRON. LE.O. OR. IRON. GT.NMERS) GO TO 130 


25 




IWRITE»2 


26 




IF (ICOL.LE.O.OR.ICOL.GT.NMEAS) GO TO 130 


27 




IHRITE-3 


28 




IF (IROW.GT. ICOL) GO TO 130 


29 




IF (IC0L.E0.3) RETURN 


30 




IHRITE»4 


31 




IF (ICOL. EO. IRON. AND. TSTORl.NE. 1.000) GO TO 130 


32 




IHRITE*5 


33 




IF (TST0R1.lt. -1.000) GO TO 130 


34 




IWRITE-6 


35 




IF (TSTORl.GT. 1.000) GO TO 130 


36 




RETURN 


37 


120 


CONTINUE 


38 




IF (IHRITE.NE. l.AND. I0IN.NE.5) WRITE (6,280) IROW 


39 




60 TO (190,140,150) , IWRITE 


40 


130 


CONTINUE 


41 




GOTO (190,200,210, 160, 170, 180) , IWRITE 


42 


140 


CONTINUE 


43 




WRITE (6,230) 


44 




60 TO 220 


45 


150 


CONTINUE 


46 




WRITE (6,240) 


47 




GO TO 220 


48 


160 


CONTINUE 


49 




WRITE (6,250) IROW, ICOL 
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GO TO 220 
170 CONTINUE 

WRITE (6.260) IROW.ICOL 
GO TO 220 
180 CONTINUE 

WRITE (6.270) IROW.ICOL 
GO TO 220 
190 CONTINUE 

WRITE (6.290) IROE.NMEflS 
GO TO 220 
200 CONTINUE 

WRITE (6.300) ICOL.NMEfiS 
GO TO 220 
210 CONTINUE 

WRITE (6.310) IROW.ICOL 
GO TO 220 
220 CONTINUE 
IFflTflL«l 
RETURN 

230 FORMAT (lOX, '«««REflS0N115O STflNOflRD OEVIflTIONS LESS THAN ZERO') 

240 FORMAT (lOX. '«nmREAS0N1200 STO. GREATER THAN MEAN VALUE.') 

250 FORMAT dOX. '*<««REAS0N1250 DIAGONAL ELEMENT IN ROW M 2. /, 1 OX. 
ICOLUMN '.12.' IS NOT EQUAL TO 1.0') 

260 FORMAT (lOX, '*<««REAS0N1300 ENTRT IN COLUMN '.12,' ROW M2./.10X.' 

less than -1.0') 

270 FORMAT (lOX. '««mREAS0N1350 ENTRT IN COLUMN '.12,' ROW M2./,10X,' 

ImmkIS greater than -^1.0') 

280 FORMAT (8X, 12) 

290 FORMAT (lOX. '^*«*«REAS0N1500 ROW M2.' IS LESS THAN ZERO OR GREATER 
1 THAN ',12) 

300 FORMAT (1 OX. '-hmREASON 1 600 COLUMN M2.' IS LESS THAN ZERO OR GREA 
ITER THAN '.12) 

310 FORMAT dOX. '« hmREAS 0N1700 ROW '.12. 'COLUMN ',12,' IS OUTSIDE OF 
lUPPER',/, lOX, ’»< m«TRIANGULAR REGION OF CORRELATION MATRIX.') 

ENO 
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SUBROUTINE SRVE 

COMMON/MCSM/CORR (78) .flNTM (2. 12) .V (25) .NMEflS.KN, 

1 SQROOT (78) ,X (25) • lU.NTMEflS. I COUNT, 

2 INflME (10) , ITMERS,NlMEflS,MftXSflM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICflRO 

REWIND NTMEflS 

WRITE (NTMEfiS) NMEfiS, KN, I NfiME 

WRITE (NTMEflS) ( (flNTM ( 1 0, JO) ,10=1,2) ,JO=l,NMEflS) 

WRITE (NTMEflS) (CORR ( 1 0) , 1 0= 1 . KN) 

REWIND NTMEflS 
WRITE (6,100) 

RETURN 

100 FORMAT!' MEflSEUREMENT DATA SAVED ON FILE ANTSTfl OR TAPE 2*) 
END 
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SUBROUTINE SlMSUB (MTTPE) 

COMMON/MCSM/CORR (78) ,RNTM (2, 12) . V (25) .NMERS.KN, 

1 SQROOT (78) ,X (25) ,1U,NTMERS. ICOUNT, 

2 INRME (10) , ITMERS,NIMERS,MRXSRM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICRRD 
DIMENSION T (45) 

1J»0 

DO 100 I=l,NMEflS 
V (I) =0.0 

CRLL RRN02 (RNORM, MTTPE. QURN, IJ) 

100 T(1)»RN0RM 

DO no I = 1.NMERS 
DO no J=1.I 
K = N0X (J.l) 

no V (1) =v(i) (j) nSqroot (K) 

RETURN 

END 
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SUBROUTINE SQRMX (NERROR) 

COMMON/MCSM/CORR (78) ,flNTM(2, 12) ,V (25) ,NMEflS,KN, 

1 SQROOT (78) ,X (25) , lU^NTMEflS, ICOUNT, 

2 I NAME (10) , ITMEflS,NlMEflS,MflXSfiM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICARD 
DIMENSION ZX (7) 

INTEGER ZX 
NERR0R»0 

DO 120 I=l,NMEflS 
ZX (M) =»I 
N»NDx (I, n 

IF (SQROOT (N) .LT.O) GO TO 130 
SO=SQRT (SQROOT (N) ) 

DO 100 J*I,NMEflS 
N»NDX (I, J) 

SQROOT (N) =SQROOT (N) /SQ 
100 CONTINUE 

IF (I.EQ.NMEflS) NERR0R*0 
IF (I.EQ.NMEflS) RETURN 
I1»I^1 
ZX(5) -11 

DO 110 I2=lKNMEflS 
DO 110 J2*I2,NMERS 
N=*N0X (12, J2) 

L-NDX (1,12) 

K«NDX (I, J2) 

SQROOT (N) =SQROOT (N) -SQROOT (L) ^SQROOT (K) 

ZX (6) =12 
ZX (7) =J2 
no CONTINUE 
ZX (1) =N 
ZX (2) =L 
ZX (3) »K 
120 CONTINUE 
130 CONTINUE 
NERR0R=1 
HRITE (6,180) 

IMO CONTINUE 

HRITE (6,190) 

CALL INPT (M,l,NVflLSZ,IflNS,XVflL,IEOF) 
1E0FP1=IE0F+1 

GOTO (170,160,150,170,160,170), lEOFPl 
150 CONTINUE 

HRITE (6,200) 

HRITE (6,210) (ZX (I) , 1 = 1,7) , SQROOT (N) 

160 CONTINUE 
RETURN 

170 CONTINUE 

HRITE (6,220) 

GO TO IMO 
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180 FORMAT (lOX, •« n«SQRMX 05I ANTHROPOMETRIC CORRELATION VALUES ARE INV 
lALlOM 

190 FORMAT (lOX,*— ENTER CARRIAGE RETURN TO CONTINUE PROCESSING \ 10 
IX,’-- ENTER ? TO OISPLAT SQRMX VARIABLES.’) 

200 FORMAT (lOX, ’« m«SQRMX 062 AN INDEFINITE QUANTITY HAS BEEN OETECTED 
IBT’,/, lOX* ’* o<kSQRMX. variable LIST WITHIN SQRMX: ’,/, lOX ’ N L 
2 K 1 11 12 J2 SQROOT ’,/, 10X,U4 (1H-) ) 

210 FORMAT (8X, 4 (3X, 12) ,4X, 3 (12, 3X) , F6. 3,/, 10X,44 (1 H-) ,//) 

220 FORMAT UOX, ’«h«SQRMX060 INVALID RESPONSE^hh ’) 

END 
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\ FUNCTION NDX (1. J) 

2 IF (l.GT.J) GO TO 100 

3 N0X=1» (JxJ-J) /2 

V RETURN 

5 100 NDX*J+ (1x1-11 /2 

6 RETURN 

7 END 
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FUNCTION UNFORM (MTTPE) 

COMMON /MCSM/ CORR (78) , flNTM (2, 12) , V (25) , NMERS, KN, 
ISOROOT (78) , X (25) , lU.NTMEflS, ICOUNT, 

21NflME (10) ,ITMEflS,NIMEflS.MflXSflM,NTEMP 
COMMON/MISC/IOIN, lOOUT, ICflRO 
DflTfl IX/ 1677721 3/, 0/281474976710655./ 

IF (MTTPE. EQ. 32) GO TO 100 

IU=1U*<IX 

UNFORM-IU/C 

RETURN 

100 1U=IU«65539 

IF (IU.6E.0) 60 TO 110 
1U*IU-^2147483647^1 
110 UNFORM-IU 

UNFORM=UNFORMk. 465661 3E-09 

RETURN 

END 
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SUBROUTINE HELP (I) 

RETURN 

END 
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99999 

100000 

100001 

100002 

100003 

looooy 

100005 

100006 

100007 

100008 

100009 

100010 

100011 



BLOCK OflTfl 

COMMON/MCSM/CORR (78) .fiNTM(2,12) . V (25) .NMEfiS.KN, 

1 SQROOT (78) .X (25) , lU.NTHERS, ICOUNT, 

2 INflME (10) . ITMERS.NlMERS.MfiXSfiM.NTEMP 
COMMON/MISC/IOIN. lOOUT. ICRRD 

DRTR 101 N, I00UT,NTEMP.NTMEflS.NlMEfiS/5.6.3,il.7/ 
DRTfi MfiX3RM/^00/ 

DRTR NTMERS.NlMEfiS/2,4/ 

1. lCOUNT/0/ 

DRTR RNTM/2M*<O./.C0fiR/78HO./.INRME/i4x‘lH / 
END 
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SUBROUTINE PUNRT 

COMMON/MCSM/CORR (78) ,QNTM (2, 12) .^(25) .NMEflS, KN, 

1 SQROOT (78) ,X (25) , lU,NTMEflS, ICOUNT, 

2 INflME (10) , ITMEflS,NlMEflS,MflXSflM,NTEKP 
COMMON/MISC/IOIN, lOOUT* ICftRD 

REWIND NTEMP 
REWIND 9 

2 CONTINUE 

REfiO (NTEMP, ENO^IMO) lNRME,NMEfiS 
WRITE (9,271) iNRME.NMEflS 

4 CONTINUE 

READ (NTEMP, END=370) (V (K) ,K*l,NMEfiS) 

WRITE (9,271) (V(K) ,K=l,NMEflS) 

271 FORMAT (20A4) 

GO TO 4 

370 CONTINUE 
ENOFILE 9 
WRITE (6,272) 

272 FORMAT (lOX, DATA SAVED ON FILE 9M 

REWIND NTEMP 

RETURN 

140 CONTINUE 
WRITE (6, 141) 

141 FORMAT (lOX, ' ERROR DATA SET D0E3 NOT EXiSTM 
RETURN 

END 
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